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Abstract 

Recent years have seen increased theoretical and experimental effort towards 
the first-ever detection of cosmic-ray antideuterons, in particular as an indirect 
signature of dark matter annihilation or decay. In contrast to indirect dark mat¬ 
ter searches using positrons, antiprotons, or y-rays, which suffer from relatively 
high and uncertain astrophysical backgrounds, searches with antideuterons ben¬ 
efit from very suppressed conventional backgrounds, offering a potential break- 
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through in unexplored phase space for dark matter. This article is based on 
the first dedicated cosmic-ray antideuteron workshop, which was held at UCLA 
in June 2014. It reviews broad classes of dark matter candidates that result in 
detectable cosmic-ray antideuteron fluxes, as well as the status and prospects 
of current experimental searches. The coalescence model of antideuteron pro¬ 
duction and the influence of antideuteron measurements at particle colliders 
are discussed. This is followed by a review of the modeling of antideuteron 
propagation through the magnetic fields, plasma currents, and molecular ma¬ 
terial of our Galaxy, the solar system, the Earth’s geomagnetic field, and the 
atmosphere. Finally, the three ongoing or planned experiments that are sensi¬ 
tive to cosmic-ray antideuterons, BESS, AMS-02, and GAPS, are detailed. As 
cosmic-ray antideuteron detection is a rare event search, multiple experiments 
with orthogonal techniques and backgrounds are essential. Therefore, the com¬ 
bination of AMS-02 and GAPS antideuteron searches is highly desirable. Many 
theoretical and experimental groups have contributed to these studies over the 
last decade, this review aims to provide the first coherent discussion of the rele¬ 
vant dark matter theories that antideuterons probe, the challenges to predictions 
and interpretations of antideuteron signals, and the experimental efforts toward 
cosmic antideuteron detection. 

Keywords: antideuteron, cosmic ray, dark matter, coalescence, propagation 


Contents 


1 Introduction |3] 

2 Prospects for dark matter detection with antideuterons |6] 

2.1 Dark matter annihilation . [7] 

2.2 Decaying dark matter. [10] 

2.3 Prospects for dark matter detection with antihelium. |Tl] 

3 Antideuteron production 1131 

3.1 The coalescence model. US! 

3.1.1 Basic concepts. 1131 

3.1.2 Determination of the coalescence momentum pq . [TS] 

3.1.3 Hadronization model dependence. |T3] 

3.1.4 Alternatives to the coalescence model .... . HH 

3.2 Gollider measurements of antideuteron production. |20| 

3.2.1 BaBab. . ED 

3.2.2 ALIGE . E3| 

3.2.3 NA61/SHINE. E3] 

3.2.4 The PANDA fixed target experiment. EH 

3.3 The coalescence model for antihelium . Ell 


2 

















4 Antideuteron propagation 1281 

4.1 Transport in the Galactic environment. [21] 

4.1.1 Cosmic-ray transport equation. [291 

4.1.2 Two-zone diffusion model. |32] 

4.1.3 Recent experimental constraints on diffusion model pa¬ 
rameters . nn 

4.1.4 Effect of dark matter density profiles. US] 

4.1.5 Galactic transport of antihelium. [36] 

4.2 Transport in the solar environment. ]36] 

4.2.1 Transport equation. HI] 

4.2.2 Force-field approximation. [H] 

4.2.3 Numerical simulation . 117] 

4.3 Geomagnetic deflection. [39] 

4.4 Atmospheric influence. ]40] 

4.5 Antideuteron interactions in particle physics detectors. [H] 

4.6 Extragalactic transport of antideuterons. SO] 

5 Experiments for the detection of cosmic ray antideuterons 1411 

5.1 BESS-Polar II: The Balloon-borne Experiment with Supercon¬ 
ducting Spectrometer . SI] 

5.1.1 BESS-Polar II instrument design. [32] 

5.1.2 BESS-Polar II antideuteron identification. [43] 

5.1.3 BESS-Polar II current status. ]34] 

5.2 AMS-02: The Alpha Magnetic Spectrometer. [33] 

5.2.1 AMS-02 instrument design . [31] 

5.2.2 AMS-02 antideuteron identification. [31] 

5.2.3 AMS-02 current status. ]3H] 

5.3 GAPS: The General Antiparticle Spectrometer. [39] 

5.3.1 GAPS instrument design . [10] 

5.3.2 GAPS antideuteron identification. [H] 

5.3.3 GAPS current status. [12] 

5.4 Experimental path forward . [H] 

6 Conclusion 1541 

1. Introduction 

The fields of cosmic-ray, antimatter, and dark matter physics have developed 
together over the last fifty years. Antideuterons, which are nuclei composed of 
one antiproton and one antineutron, were identified in the early 1960s among 
the secondaries from nuclear targets exposed to the proton beams of the CERN 
Proton Synchrotron [T] and Brookhaven Alternating Gradient Synchrotron |2]. 
By this time deuterons, which are nuclei composed of one proton and one neu¬ 
tron, had already been identified in nuclear emulsions exposed to cosmic rays jS] . 
The production of deuterons and antideuterons in collider experiments, however, 
allowed for tests of models describing their formation [1111[S1[7]. From this. 
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the picture of a two-stage formation process emerged, with the production of a 
nuclear cascade within which several nucleons interact to form a nucleus. 

By the 1970s, the existence of dark matter had been established on galactic 
scales by precision galaxy rotation measurements pin]. The presence of enough 
hidden, heavy, normal objects to account for this dark matter has now been 
ruled out by astronomical surveys and cosmological simulations m. which con¬ 
firm that dark matter is roughly five times more abundant in the universe than 
baryonic matter m- Current experiments instead focus on searching for en¬ 
tirely new particles that will make up this mass. Among the most theoretically 
well-motivated candidates for such a particle is the weakly interacting massive 
particle, or WIMP. 

It is generally acknowledged that the identification of dark matter will re¬ 
quire the interplay of direct searches (which probe scattering cross sections 
through the recoil of dark matter off of target nuclei), indirect searches (which 
probe annihilation cross sections and decay lifetimes through the detection of 
annihilation or decay products), and collider searches (which probe production 
cross sections through the products of high-energy particle collisions). There 
are currently about twenty operating or planned direct search experiments 
[El [HUH]. The current best exclusion limits come from the LUX and CDMS ex¬ 
periments [HiiiKn]. However, these limits conflict with earlier CDMS-II [T5] . 
DAMA/LIBRA [H HOI I1I1I21 US] , and CoGeNT HHSJ [H HZl US] results, 
which are consistent with light (« 10 GeV) dark matter. The results of CDMS- 
II and CoGeNT can be made compatible with the limits of LUX [TS] , CDMSlite 
[TS] and XENON |29] if dark matter is not isospin-invariant [SO] jHI] [32]. In 
any case, the high backgrounds and low signal amplitudes of direct detection 
experiments make analysis, particularly at low WIMP mass, challenging, and 
direct experiments optimized for low-energy recoils, such as PIGO |33) . are still 
orders of magnitude less sensitive for the spin-independent case than LUX or 
GDMS. 

If dark matter was in thermal equilibrium in the early universe, and froze 
out when the temperature dropped due to expansion, it is a natural assumption 
that dark matter particles are able to interact with each other and produce 
Standard Model particles. Indirect searches exploit possible kinematic differ¬ 
ences between the production of cosmic rays through dark matter and standard 
astrophysical processes to identify dark matter signals. Gosmic-ray antiparti¬ 
cles are ideal candidates for indirect searches, but are often challenging due 
to high or uncertain levels of astrophysical background. Recently, results from 
PAMELA |Mj, Fermi [3^, and AMS-02 jSS] [ST] [35] show a structure in the 
positron data that might be interpreted as originating from dark matter with 
TeV-scale mass. Antiproton data do not show such a significant feature [3^14(1] . 
These null antiproton results instead have been used to constrain a variety of 
dark matter models 1111111113], though their statistical accuracy at low energies 
is limited. 

Antideuterons may be generated in dark matter annihilations or decays, 
offering a potential breakthrough in unexplored phase space for dark matter. 
In this process, annihilation or decay proceeds via particle physics, produc- 
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ing particles that then hadronize to an antiproton, p, and antineutron, h, pair 
that coalesce to form an antideuteron, d. The unique strength of a search for 
low-energy antideuterons lies in the ultra-low astrophysical background for this 
channel and was pointed out for the first time by the authors of [Hj. By now, 
the topic has been further discussed in more depth for more than a decade, 
e.g., in [11113 SSI EZl HE]- The dominant conventional sources for secondary 
(background) antideuteron production are cosmic-ray protons or antiprotons 
interacting with the interstellar medium |47j . However, the high threshold en¬ 
ergy for antideuteron production and the steep energy spectrum of cosmic rays 
mean there are fewer particles with sufficient energy to produce secondary an¬ 
tideuterons, and those that are produced have relatively large kinetic energy. 
In addition, the low binding energy of antideuterons makes energy loss through 
collisions difficult. 

Fig.[l]shows the antideuteron flux expected from three benchmark dark mat¬ 
ter scenarios. These dark matter candidates include a lightest supersymmetric 
particle (LSP) neutralino from the minimal supersymmetric model (MSSM), a 
5D warped GUT Dirac neutrino (LZP), and an LSP gravitino. The expected 
secondary/tertiary background |49j is also shown. This figure reveals why low- 
energy antideuterons are such an important approach: the flux from a wide 
range of viable dark matter models exceeds the background flux by more than 
two orders of magnitude in the energy range below 0.25 GeV /n, and by more 
than an order of magnitude up to 1 GeV/n. 

The large signal-to-background ratio for low-energy antideuterons does not 
rely on any boosting mechanisms due to, e.g., dark matter dumpiness, Sommer- 
feld enhancement, or large Galactic Halo size. Introducing such effects would 
further increase the signal strength. This is in contrast to dark matter signal 
predictions for positrons, antiprotons, and y-rays, which typically constitute a 
small component on top of the astrophysical background. Sec. [previews various 
dark matter and other beyond the Standard Model candidates that produce 
an antideuteron signal evading current experimental bounds, but within the 
sensitivity of ongoing or planned cosmic-ray antideuteron experiments. 

The coalescence model describes the process through which an antideuteron 
is formed through the merging of the antiproton and antineutron. As outlined 
in Sec. several implementations of this model have been developed. The 
first analyses |44l HU |46l EH [Ml EZ] were performed within a purely analytical 
framework, while more recent ones jUj |3H1 SSI ESI ESI ESI ESI 1^ El] exploit an 
event-by-event Monte Garlo approach. These coalescence models are then tuned 
to collider measurements of deuteron and antideuteron production, as described 
in Sec. 13.21 

The predicted cosmic antideuteron flux arriving at the Earth relies on the 
modeling of charged-particle transport both in the Galactic medium and in the 
heliosphere. These transport models are discussed in Sec. [^ The modeling 
of Galactic transport is the main source of theoretical uncertainty, while solar 
modulation is relevant to shaping the low-energy tail of the predicted flux. 

Despite the large signal-to-background ratio, any search for cosmic-ray an¬ 
tideuterons is a rare event search, and multiple experiments with complementary 
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Figure 1: Predicted antideuteron flux as a function of kinetic energy per nucleon for a 30 GeV 
neutralino, a 40 GeV extra-dimensional Kaluza-Klein neutrino, and a 50 GeV gravitino |45l 
14611501149| . The antideuteron limits from BESS are shown ED. along with the projected 
sensitivities of AMS-02 for the superconducting-magnet configuration m after 5 years of 
operation and GAPS aft er th ree 35-day flights [531 I54| . The MED Galactic propagation 
scenario is assumed (Sec. [4T) |. These predictions use a coalescence momentum that is set 
to 195 MeV (Sec. and the Einasto dark matter density profile (Sec. [47l.4[ |. For the solar 
modulation parameters see Sec. |4.2| 


techniques and backgrounds will be necessary to build confidence in any detec¬ 
tion. Cosmic antideuterons have not yet been detected, with the current best 
flux upper limits provided by the BESS experiment m- More sensitive limits 
will be provided by the AMS-02 experiment [531 IMl EH [55] , which is currently 
taking data onboard of the International Space Station. The General Antiparti¬ 
cle Spectrometer (GAPS) experiment [571155115511751 ITT] . which is proposed for 
several Antarctic balloon campaigns, will provide essential complementary sen¬ 
sitivity to both antideuterons and antiprotons in an unprecedented low-energy 
range. Sec. [^presents the experimental sensitivities and status of the BESS, 
AMS-02, and GAPS antideuteron searches. 

A number of reviews on the status of dark matter searches have been pub¬ 
lished in recent years, with discussions related to antideuteron searches appear¬ 
ing in, e.g., the following [72l|73l|74l|75]. 

2. Prospects for dark matter detection with antideuterons 

Many dark matter models are capable of producing an antideuteron flux that 
is within the reach of currently operating or planned experiments. Sec. |2.1| dis- 
cusses antideuterons resulting from dark matter self-annihilations, while Sec.|2.2| 
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discusses antideuterons from decays of long-lived dark matter. The purpose of 
the discussion is not to detail every theoretical model that can result in a de¬ 
tectable flux of antideuterons; instead, the following review demonstrates that 
measurements of cosmic-ray antideuterons are sensitive to a wide range of theo¬ 
retical models, probing dark matter masses from 0(1 GeV) to 0(1 TeV). Anti¬ 
helium nuclei can be produced by many of the same beyond the Standard Model 
processes as antideuterons, providing an even higher signal-to-background ratio 
at low energies, but also significantly lower predicted flux levels. The prospects 
and challenges of future dark matter searches using antihelium signatures are 
discussed in Sec. 12.31 

Before starting the discussion, it is important to point out the dominant 
uncertainties on the antideuteron flux predictions presented here, which will 
be more deeply addressed in Sec. and These uncertainties are due to the 
hadronization and coalescence models used to describe antideuteron formation, 
as well as the propagation models used to describe antideuteron transport in 
the Galactic and solar environments. The antideuteron fluxes presented in this 
section assume a conservative boost factor, due to dark matter clumps in the 
Galactic halo, of / = 1. However, a boost factor of / = 2-3, as is consistent 
with current theoretical expectations, would increase dark matter fluxes by a 
factor / over those discussed below [JSl [TSl [77]. Such a boost factor is only 
relevant for dark matter annihilation, as dark matter decay depends linearly on 
the dark matter density. 

It is also vital to note that every process capable of producing antideuterons 
will also produce a much larger flux of antiprotons. Thus, the strength of 
any prospective antideuteron signature from dark matter is constrained by the 
(current) non-observation of a signature in the cosmic antiproton spectrum. 
However, the ultra-low background for low-energy antideuteron searches still 
opens new phase space for dark matter observations. Furthermore, detecting 
deviations from the astrophysical antiproton flux requires very high statistics. 
Antideuterons provide an additional search channel with very strongly sup¬ 
pressed astrophysical backgrounds compared to antiprotons, and thus can act 
as an essential probe to confirm or rule out potential deviations in the antiproton 
spectrum due to processes like dark matter annihilation or decay. On the other 
hand, a non-detection of a signal above background in the antiproton channel 
might also just be a consequence of experimental limitations and can serve as 
additional motivation to pursue antideuteron searches. 

2.1. Dark matter annihilation 

It is common to illustrate experimental sensitivity to annihilating dark mat¬ 
ter by first assuming an annihilation cross section into a particular Standard 
Model channel, without considering the details of the underlying dark matter 
model. Analytical models or Monte Carlo simulations can then be used to cal¬ 
culate the resulting flux level of different cosmic-ray species. Since larger dark 
matter particle mass implies lower dark matter particle density, the lowest anni¬ 
hilation cross section that a measurement is sensitive to then scales quadratically 
with the dark matter mass. 
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Figure 2: The sensitivity of GAPS to a three standard deviation detection in the antideuteron 
channel of a dark matter WIMP annihilating into uu, shown as a function of annihilation 
cross section, (crv), and dark matter mass, muM- The red, black, and green lines indicate the 
GAPS sensitivity in the MIN, MED, and MAX Galactic propagation scenarios (Sec. [4^. The 
dark matter density profile, solar modulation parameters, and coalescence momentum are as 
defined in Fig.j^ Source: Adapted from 1421 . 


The prospects for a three standard deviation detection with the GAPS ex¬ 
periment of a dark matter WIMP annihilating into mm, as a function of anni¬ 
hilation cross section and dark matter mass, are shown in Fig. The three 
solid lines refer to the GAPS sensitivity assuming the MIN, MED and MAX 
Galactic propagation models (Sec. 4.11, while the dashed line represents the 
thermally averaged cross section for a thermal-relic WIMP. The theoretical and 
experimental uncertainty for the propagation model and coalescence scheme are 
both on the order of 0(10). The experimental sensitivity reaches well below 
the thermal-relic cross section for dark matter masses below todm = 20 (100) 
GeV for the MED (MAX) propagation scenario. It is important to stress that 
recent measurements disfavor the MIN model, as will be further discussed in 
Sec. |4.1.3[ Antideuteron searches thus both provide complementary sensitiv¬ 
ity to direct search experiments, which are most sensitive to intermediate-mass 
dark matter, and provide increased sensitivity to light dark matter models, 
where direct search techniques are challenging. For example, light candidates 
in the next-to-minimal supersymmetric model (NMSSM) still evade direct, col¬ 
lider, and Fermi dwarf galaxy bounds [75] , but would be probed by ongoing and 
proposed antideuteron experiments. An example of the sensitivity to standard 
supersymmetric particles is shown in Fig.[^ which presents the signal predicted 
for a 30 GeV supersymmetric neutralino annihilating into bb. 

This signal is within the detectable range of both GAPS and AMS-02, and is 







motivated by the diffuse y-ray excess observed surrounding the Galactic Center 
by Fermi m- The excess shows a high significance and the interpretation is 
ongoing [80l [U US [83l HH |85] . 

In addition, the authors of [53] state that the gluon-gluon decay channel 
of dark matter could produce an antideuteron flux in the detectable range for 
masses below 0(100 GeV). Also shown in Fig. [^is the case of a right-handed 
Kaluza-Klein neutrino of warped 5-dimensional grand unified theories with a 
conserved Z 3 parity (LZP) |45l|86l[87l[88] , which for a mass of 40 GeV has a ma¬ 
jor annihilation channel into Z bosons. The authors of [45] interpret especially 
the LZP mass range close to the Z pole of 40-50 GeV as an interesting bench¬ 
mark case for WIMP pair annihilation into Z that still produces the correct 
thermal relic density. 

Antideuteron signals from annihilation of heavy (tudm = 0.5-20 TeV) dark 
matter may also be detectable by AMS-02 and GAPS, for the MAX propagation 
scenario [56] . These multi-TeV mass dark matter particles are motivated by the 
possible high-energy cosmic positron excess observed by AMS-02 [551 133 I3H] • 
Accounting for the lack of a similar excess in the available antiproton data, it 
requires both high dark matter mass and an enhanced annihilation cross sec¬ 
tion, such as provided by the Sommerfeld mechanism [89] . The prediction for 
such multi-TeV mass particles annihilating into bb is shown in Fig. Anni¬ 
hilation into W'^W~ would be disfavored in this scenario. Another model of 
heavy supersymmetric dark matter, also shown in Fig. [^ is pure-Wino dark 
matter m- The thermal relic dark matter density is reached for Wino dark 
matter at masses in the TeV range, which is close to the detectable range in the 
low-energy region of GAPS and AMS-02, and in the detectable range for the 
high-energy region of AMS-02. It is interesting to note that the antideuteron 
spectra from quark and heavy gauge boson annihilation channels are qualita¬ 
tively different. Under the assumption that dark matter is moving slowly and 
that the dark matter halo and the Earth’s rest frame are roughly the same, the 
hadronization in gauge boson channels is happening in a boosted frame while 
the hadronization in quark channels is essentially happening in the Earth’s rest 
frame. As a result, the shapes of the antideuteron spectra originating from 
quark annihilation channels only vary slightly with mass and the gauge boson 
channels show a stronger dependence. 

Recent Higgs portal dark matter models discuss the possibility of dark mat¬ 
ter annihilating via Higgs exchange processes [3T]I35]I351[M]. In these models, 
the scalar field S plays the role of cold dark matter in the universe with a 
scalar gauge singlet mass, ms- There are three mass regions where it is possible 
to reproduce the correct dark matter relic abundance. The low-mass region 
(ms ^ 50GeV), in which the dominant annihilation channel is into bb pairs, is 
completely ruled out by direct detection experiments and the bound on the in¬ 
visible decay width of the Higgs [35]. The resonant region (ms « 50 — 70 GeV), 
in which dark matter annihilation proceeds via the on-shell exchange of the 
Higgs boson, is severely constrained by cosmic antiproton [96] and y-ray m 
measurements. The high-mass region (ms ^ 70 GeV), in which the domi¬ 
nant annihilation channels are into electroweak gauge bosons, two Higgses (if 
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Figure 3: Predicted antideuteron flux for annihilation of dark matter with = 5, 10, 

20 TeV |56| (blue lines, top to bottom) into bb with enhanced annihilation cross sections 
((crD}5TeV = 3- 10“^^cm3/s, ((T'u}lOTeV = 7- lO^^^cm^/s, (o-'u}20TeV = 20- 10“^^ cm®/s). 
The predicted antideuteron flux from pure-Wino dark matter m (solid green line) with 
= 0.5 TeV, (uv) = 4.82 • 10~^^cm^/s with Sommerfeld enhanced annihilations into 
W'^W~. The MAX propagation model is used for all predictions. 


ms > 2mh « 252 GeV), and ti pairs (if ms > 2mt « 347 GeV), are still only 
marginally disfavored by direct and indirect detection experiments [96], and 
for ms > 100 GeV one can obtain a good dark matter candidate. This region 
extends up to the TeV scale. Dark matter candidates of this type with inter¬ 
mediate dark matter masses are able to produce an antideuteron signal well 
above the astrophysical background, but not yet within the reach of the current 
experiments. 

2.2. Decaying dark matter 

Dark matter may also exist in the form of very long-lived, but unstable 
particles whose decays can produce antideuterons |98] . A particularly well- 
motivated and well-studied case of decaying dark matter is the gravitino in 
supergravity models with i?-parity violation |99l IIOOL llOlj . In the most general 
case, i?-parity violation adds one bilinear and three trilinear renormalizable 
operators to the superpotential [102] : 

^^ XijkLiLjEk + + X'k^UiDjDk , ( 1 ) 

where Li and Qi are the up-type Higgs, lepton and quark SU(2) doublet 
superfields, respectively, and Ei^ Ui and Di are the lepton, up-type quark and 
down-type quark SU(2) singlet superfields, respectively. Summation over gen¬ 
eration indices i,j,k and suppressed gauge indices is assumed. The first three 
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operators violate lepton number and the last operator violates baryon number. 
Moreover, there could be additional bilinear and trilinear i?-parity violating 
operators in the soft supersymmetry breaking terms. Proton stability is not 
violated as long as only specific subsets of those operators are present in the 
theory m- A theory with bilinear ii-parity violation only violates lepton 
number conservation, while theories with trilinear i?-parity violation usually 
assume the presence of either only lepton number violating operators or only 
baryon number violation. Since the gravitino decay rate is quadratically sup¬ 
pressed by the Planck mass, 1/Mpj, moderately small i?-parity violating param¬ 
eters are sufficient to make the gravitino long-lived on cosmological time scales 
(r 3> 4 X 10^"^ s). 

Taking into account constraints on the gravitino lifetime derived from the 
PAMELA antiproton observations m , i.e. T > 10^®-10^®s [43], antideuteron 
signals from gravitino decays into the channels ■jVi, Zvi, W(.i, and hvi via 
bilinear R-parity violation |105l llObj exceed the astrophysical background at low 
energies, but are not within the detectable region of AMS-02 or GAPS [sniisi 
snz]. Gravitino decays via trilinear i?-parity violating operators, however, may 
be detectable jsnmnHj- In these models, the gravitino typically decays into three 
final state particles: Vidjdk, iiUjdk, Uidjdk [109j . The purely leptonic operator 
LiLjEk leads to final states Vi£jik and thus does not produce any antideuterons. 
A particularly relevant case is the baryon number violating operator UiDjDk 
that triggers gravitino decay into three quarks. The production of three quarks 
in the hard process leads to a larger number of antiprotons and antineutrons 
in the final state. While this leads to stronger constraints from antiproton 
observations, the increased probability for antideuteron formation leads to an 
enhancement of the expected antideuteron flux compared to the antiproton flux. 
Therefore, this particular trilinear i?-parity violating gravitino decay channel 
allows for a flux that is within the reach of AMS-02 and GAPS while obeying 
the antiproton constraints m- This is demonstrated for the case of a 50 GeV 
gravitino in Fig. 

Note that since the antideuteron flux only drops linearly with increasing 
gravitino mass (as opposed to for annihilating dark matter), the grav¬ 

itino decay signal may dominate over the background of secondary antideuterons 
at high energies as well. This opens the possibility of using high-energy an¬ 
tideuterons (i.e. T > 100 GeV) to look for signals from heavier decaying dark 
matter. 

2.3. Prospects for dark matter detection with antihelium 

Due to the larger energy threshold necessary to produce secondary anti¬ 
helium nuclei (31 Wp = 29.09 GeV, instead of 17 mp = 15.95 GeV to produce 
antideuterons), antihelium signatures of dark matter annihilation or decay have 
an even higher signal-to-background ratio of 10^ — 10^ (Fig. at low energies 
than antideuteron signatures (about 10^) |11111112] . However, in these mod¬ 
els also the antihelium formation from dark matter is suppressed compared to 
antideuteron formation. 
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Figure 4: Predicted antihelium flux from various dark matter models. The red regions show 
the flux from of a 40 GeV dark matter particle annihilating into bb nn] with (crv) = 3 • 
10“^®cm^s~^. The solid red band corresponds to coalescence momentum pQ = 195MeV 
and the hatched red band to po = 300 MeV (Sec. Only those configurations that are 
compatible with PAMELA antiproton constraints are considered. The predicted antihelium 
flux of a 100 GeV dark matter particle annihilating into W'^W~ with ( 0 - 1 ;} = 3-10 cm^ s ^ 
is shown in blue For both cases, the band width reflects variations due to different 

propagation model parameters. The astrophysical background is represented by the green 
solid line. The projected sensitivity of AMS-02 after 5 years of operation is shown by the light 
green region. The gray regions correspond to the bounds imposed by previous BESS [m], 
PAMELA [115| . and AMS-01 [116| measurements (the helium flux from PAMELA was used 
to translate from the He/He results given in these references). The dark matter density profile 
is Einasto. 


Fig.i shows the predicted fluxes for a 40 GeV dark matter particle anni¬ 
hilating into bb and a 100 GeV dark matter particle annihilating into W~^W~. 
The thermal-relic annihilation cross section is used for both predictions. The 
width of the shaded bands is defined by the variation in flux predicted using 
a range of Galactic propagation model parameters that produce an antipro¬ 
ton flux compatible with PAMELA measurements. As discussed in Sec. |3.3[ 
there is a large uncertainty due to the choice of the coalescence momentum, 
Po- This is illustrated by the two red bands corresponding to po = 195 MeV 
and Po = 300 MeV for the bb annihilation channel. The classical coalescence 
model HIT] predicts a flux dependence proportional to p® for antihelium, as 
opposed to Po for antideuterons. Thus constraining the value of the coalescence 
momentum is critical to constraining antihelium flux predictions. The sensitiv¬ 
ity of AMS-02 with 5 years of integration time m is indicated by the light 
green region. These antihelium signatures are outside the current sensitivity of 
both AMS-02 and GAPS, but motivate future, more sensitive experiments. 
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3. Antideuteron production 


Although antideuterons were discovered for the first time in 1965 mm, 
their formation is not well understood. However, understanding antideuteron 
production is crucial for the interpretation of the cosmic-ray data, which im¬ 
pacts both the antideuteron background expectation from interactions of pri¬ 
mary/secondary cosmic rays with the interstellar medium as well as the forma¬ 
tion in the aftermath of dark matter annihilations or decays. In recent years 
event-by-event Monte Carlo techniques have been used to show that the an¬ 
tideuteron flux prediction carry an uncertainty of about one order of magnitude 
for the formation mechanism. It is currently an open question if the antideuteron 
production depends on the exact underlying process and on the available center- 
of-mass energy or if Monte Carlo generators need further refinement. This 
section will explain the theoretical status as well as future outlook of the under¬ 
standing of the formation of antideuterons (Sec. 3.1). This discussion is followed 
by a description of the recent experimental results of antideuteron production in 
hadronic decays of the T resonance from BaBar (Sec. 3.2. 1[ ) and otp-p, p-Pb, and 
Pb-Pb collisions in the ALICE experiments (Sec. 3.2.2| ). In addition, the out¬ 
looks for the operational fixed target expe riment NA61/SHINE (Sec. 3.2.31 and 
the future PANDA experiment (Sec. 3.2.4) are discussed. The section concludes 
by outlining the procedure for the antihelium coalescence model (Sec. 3.3). 


3.1. The coalescence model 

3.1.1. Basic concepts 

The fusion of an antiproton and an antineutron into an antideuteron is de¬ 
scribed by the coalescence model. The coalescence model is based on the sim¬ 
plifying assumption that any (anti)nucleons within a sphere of radius po in 
momentum space will coalesce to produce an (anti)nucleus. For the case of 
deuterons, the condition means that any pn-pair with Ap < po will coalesce 
to produce a deuteron. The coalescence momentum po is a phenomenological 
quantity, and has to be determined through hts to experimental data. The 
model was first proposed by Schwarzschild and Zupancic [7] for the production 
of deuterons, tritons and helium nuclei in fixed target scattering experiments. 
This has since been adopted for the description of antideuteron production in 
collider experiments and in the context of indirect dark matter detection. 

One of the strengths of this simple model is that under the assumption of 
isotropic and uncorrelated proton and neutron momentum distributions, ana¬ 
lytic expressions can be found for the spectra of different nuclei, given in terms 
of the nucleon spectra and the coalescence parameter pq. Before the onset of 
Monte Carlo event generators, such analytic expressions were also entirely neces¬ 
sary for application of the coalescence prescription. The antideuteron spectrum 
can in this approach be given in terms of the antiproton and antineutron spectra 
through 

^ Po 1 d^ d^ 

dTj 6 mnUip Irpj _|_ 2mjTj ' 
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where rrii, Ti and dNi/dTi are, respectively, the mass, kinetic energy and differ¬ 
ential yield per event of particle i. 

In the literature, some variant of this expression will be found in most works 
before 2009, but often with a differing constant prefactor - typically a factor of 
8 - which can be absorbed into a redefinition of po- While the assumptions of 
isotropic and uncorrelated nucleon spectra might be a good approximation in 
low-energy or minimum bias nuclear interactions for which the model was made, 
these assumptions have been found not to hold in relevant elementary particle 
interactions |58j , such as dark matter annihilations or decay, and p-p collisions 
at low center-of-mass energies. This is due to the nucleons typically being 
produced in geometrically restricted QCD jets, leading to strong correlations 
between nucleons on a per-event basis. 

In order to take these correlations into account, the state of the art has since 
2009 been to apply the coalescence condition to ph-pairs on a per-event basis in 
Monte Carlo events. Unfortunately, this requires some four orders of magnitude 
more events than the isotropic approximation in order to achieve the same level 
of statistical error on predictions. As extensively discussed below, different 
event generators yield different values of po when compared to a particular 
experiment, indicating a substantial systematic uncertainty in the coalescence 
prediction. Even with the same event generator, experiments differing in energy 
and interaction type will typically give inconsistent values of Pq. 

Since momentum is not a Lorentz invariant quantity, the frame in which the 
coalescence condition is applied could also potentially be of importance. In the 
original coalescence model, the coalescence condition was applied in the rest 
frame of one of the coalescing particles, while in recent works, employing per- 
event coalescence, the condition is either applied in the center-of-mass frame of 
the antinucleons, or applied as a condition on the invariant momentum, 

-A^<Po (3) 

where = kp — k^. These two approaches turn out to be equivalent in the 
limit TTin = mp, and are in practice easily interchangeable. 

As pointed out in |48j . spatial separation should also be taken into account 
in the formation of antideuterons. Nuclear interactions typically take place on 
scales of a few femtometers. Antinucleons originating from weakly decaying 
particles, whose long lifetimes lead to macroscopic decay lengths, will hence 
be produced too far from the primary interaction vertex to have any chance 
of interacting with antinucleons from the primary collision. To take this into 
account, one can consider weakly decaying particles stable in this context. This 
considerably reduces the number of antinucleons available, and correspondingly 
increases the value of po required to reproduce experimental data. The reader 
should be aware that not all recent work takes spatial separation into account, 
and the cut on lifetime varies when used, e.g., Ref. [48] uses Imm/c while [50] 
uses lA/c. While the differences that originate from differing values of the 
lifetime cut are believed to be fairly small, one should keep this in mind when 
comparing values of po from different sources. 
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Figure 5: Results of fitting the coalescence momentum pQ to different datasets on antideuteron 
production, based on im and [50]. 

3.1.2. Determination of the coalescence momentum 

The predicted antideuteron fluxes from dark matter annihilations or decays, 
as well as from cosmic-ray spallations depend strongly on the only free param¬ 
eter of the coalescence model: the coalescence momentum pq. Currently, this 
parameter cannot be calculated from first principles, and should therefore be 
determined from fitting the predictions of the event-by-event coalescence model 
described in the previous section to the experimental results on antideuteron 
production. 

The ALEPH [m], CLEO [Iig, CERN ISR [11111110], ZEUS [HI], AL¬ 
ICE |122j . and BaBar [I23j experiments have measured antideuteron produc¬ 
tion, where the recent results of the latter two are discussed in more detail in 
Sec. 13.2.1] and The best-fit values of the coalescence momentum inferred 

from each experiment using PYTHIA |124j or HERWIG++ [125] . including the 1, 
2 and 3 standard deviations error bands, are summarized in Fig. [^ As ap¬ 
parent from the plot, when inferring the coalescence momentum using PYTHIA 
with default settings, there is no value of po that can simultaneously fit all 
the data [48]. On the contrary, the best-fit value of the coalescence momentum 
shows a dependence on the underlying process and on the center-of-mass energy 
involved. Interestingly when using HERWIG++ and simultaneously tuning po and 
three hadronization parameters to antideuteron data from ALEPH, CLEO and 
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Figure 6: Ratio between the antideuteron spectra from HERWIG++ and PYTHIA 6.4 as function 
of the scaled kinetic energy x = T/m-£,yi. Source: Figure taken from m- 


ZEUS as well as antiproton data, these three datasets can be brought in agree¬ 
ment with a single value of po [50]. However, the data from CERN ISR and 
BaBar still give rise to a significantly different pq. The hadronization model 
dependence is discussed more extensively in Sec. |3.1.3| 

As a consequence, it is currently inconclusive which value of po and which 
Monte Carlo generator is most suitable for calculating the antideuteron flux from 
dark matter annihilations or decays as well as for the production in cosmic-ray 
spallation. This uncertainty has dramatic implications for the search for cosmic 
antideuterons, due to the strong dependence of the antideuteron yield on the 
coalescence momentum, Ng cx p^. This proportionality is exact in the isotropic 
coalescence model, as can be seen from Eq. ([^, but also holds well for per-event 
coalescence for kinetic energies above po. 

3.1.3. Hadronization model dependence 

Calibration against experimental data typically yields best-fit po-values that 
are smaller than or comparable to Aqcd “ the scale at which perturbative QCD 
breaks down. The coalescence model is therefore sensitive to non-perturbative 
effects in the hadronization models of Monte Carlo event generators. This natu¬ 
rally includes (anti)baryon production rates from partons. However, due to the 
small value of the coalescence momentum p^, coalescence is also highly sensitive 
to two-particle correlations between the (anti)nucleons. 

This issue was first raised in EH, where the uncertainty from hadroniza¬ 
tion dependence was estimated by comparing antideuteron spectra produced 
in the Monte Carlo event generators HERWIG++ and PYTHIA 6.4. While PYTHIA 
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employs the Lund string hadronization model, HERWIG++ uses cluster hadroniza- 
tion. Neither PYTHIA nor HERWIG++ are properly tuned for p-p collisions at y/s « 
10 GeV, which is the dominant production channel of secondary antideuterons. 
For that purpose, the better suited Monte Carlo generator DPMJET-III |126) 
has been employed to calculate the antideuteron background flux using the 
event-by-event coalescence model [33] ■ In the comparison, the two Monte Carlos 
were calibrated separately against ALEPH data on antideuteron production in 
hadronic Z-decays m , yielding different best-fit values for the coalescence mo¬ 
mentum pq. The resulting antideuteron spectra from dark matter annihilations 
obtained with HERWIG++ and PYTHIA were found to differ by a factor of about 
2-4 for most energies, depending on the process in question, but rapidly in¬ 
creasing at high and low energies. The ratio between the HERWIG++ and PYTHIA 
antideuteron spectra as a function of the scaled kinetic energy x = T/ttidm, 
where todm is the dark matter mass, can be seen in Fig. For the XX 
annihilation channel, the point at low energies below which the uncertainty rises 
steeply was found to lie at a kinetic energy of T « 10“^ GeV, independent of 
the dark matter mass. This is below the experimentally relevant energies. For 
the XX W~^W~ channel, however, the point was observed to lie at a constant 
value of the scaled kinetic energy, x « 10“^, which means that experimentally 
relevant kinetic energies can lie in the high uncertainty region for dark matter 
masses above a few TeV. The underlying reason for this fundamental difference 
is that while the quarks, being colored (and color-connected) objects, radiate 
off increasing amounts of QCD radiation with increasing energies, the gauge 
bosons are treated as on-shell. As long as the two Ws are sufficiently boosted 
to disallow interactions between their decay products, the antideuteron spec¬ 
trum equals the sum of their spectra when decaying at rest, boosted according 
to the speed of the Ws. 

It is important to note that most of the uncertainties affecting the determi¬ 
nation of the antideuteron signal from dark matter annihilations or decays are 
strongly correlated with the corresponding uncertainties in the antiproton flux, 
which has been measured rather precisely by PAMELA |104j . For example, un¬ 
certainties stemming from unknown propagation parameters do not significantly 
affect the antideuteron-to-antiproton ratio expected from dark matter annihi¬ 
lations or decays. However, the uncertainty stemming from the two-particle 
correlations in the different hadronization models is not correlated with the an¬ 
tiproton flux, and hence deserves special attention when evaluating prospects 
for antideuteron detection. 

As the underlying mechanisms of the hadronization models are very dif¬ 
ferent, one should not expect the two-particle correlations from the two main 
hadronization models to agree. There is also no a priori reason to expect the 
two-particle correlations from one Monte Carlo to be more reliable than the 
other. However, as pointed out by m, better agreement with data might be 
achieved for both hadronization models by tuning the phenomenological param¬ 
eters in the models specifically for antideuteron production. 

Each model has a large number of free parameters, which have been tuned 
against various experimental data, typically with a focus on the ability to re- 
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Table 1: Results from a hadronization parameter fit compared to default values in HERWIG++. 
po-values are in units of MeV. 


Parameter 

Default value 

Value at xLn 

Uncertainty 

Po 

- 

143.2 

-5.5 

ClMaxLight 

3.25 

3.03 

+0.18 

-0.15 

PSplitLight 

1.20 

1.31 

+0.19 

-0.32 

PwtDIquark 

0.49 

0.48 

+0.15 

-0.04 


produce particle multiplicities and jet physics results. These general tunes do 
not take two-particle correlations into account to any significant degree, and 
(anti)proton spectra have no higher weights than other particles in the tuning. 
For the purpose of antideuteron production, the spectra of particles that do not 
contribute to the (anti)nucleon or antideuteron spectra are of no importance. A 
tuning of the hadronization parameters that only takes the relevant spectra into 
account is more likely to be able to reliably reproduce experimentally observed 
antideuteron spectra. This has previously been shown to improve the physics 
description when done with other processes, such as Higgs production |127j . A 
further advantage of this procedure is that a combined fit of po together with a 
set of hadronization parameters also allows for an estimation of the uncertainty 
on the antideuteron spectrum from the hadronization and coalescence tuning. 

In [SO], the coalescence momentum po was tuned along with three of the most 
relevant hadronization parameters in HERWIG++ against antideuteron data from 
ALEPH, CLEO and ZEUS, as well as proton/antiproton data from ALEPH and 
OPAL. The inclusion of (anti)proton data in the tuning is necessary for a con¬ 
sistent result, as one cannot expect a reliable reproduction of the antideuteron 
spectrum if the spectra of its constituents are wrong. The best-fit point of the 
tune can be found in Table and is in reasonably good agreement with the 
default hadronization parameters for HERWIG++. Eor a short description of the 
role of the various parameters, please see the original article |50|. Note that 
the parameter errors are strongly correlated, and that the given numerical val¬ 
ues therefore cannot be naively used to estimate the goodness-of-fit of other 
parameter points. 

The Po parameter was intended to describe the maximum difference in mo¬ 
menta for which two (anti)nuclei will form ions. The definition in Eq. [fallows 
Po to be determined from the measurement of (anti)nuclei spectra. However, 
Po not only describes the required difference in momenta of the coalescence 
partners, but also parameterizes a number of other effects. Eor instance, the 
available phase space for ion production depends on the available energy in the 
formation interaction. Comparing the deuteron and antideuteron production 
in p-p interactions at the same energy it is clear that a pn-pair can already be 
produced by pp —>■ pmr~^ while a pfi-pair requires at least six final-state nucleons 
(pp —>■ pnnppp). Therefore, if the coalescence approach is applied to describe 
the total deuteron and antideuteron production cross sections, po describes not 
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only the probability for merging an (anti)proton and an (anti)neutron, but also 
the probability for (anti)protons and (anti)neutrons to be produced at all. This 
is especially important for (anti)deuteron production close to the production 
threshold energy. In a similar way, the production close to the threshold also 
favors an anti-correlation of (anti)protons and (anti)neutrons for kinematical 
reasons and causes an additional phase-space suppression. The astrophysical 
background production of antideuterons in interactions of cosmic rays with the 
interstellar medium is most likely dominated by the production at the thresh¬ 
old, and thus strongly phase-space suppressed. However, the energy available 
in dark matter annihilations or decays can be much higher, and therefore the 
antideuteron production can be much less suppressed. Hence, using different 
coalescence momentum values for different dark matter masses as well as for 
different contributing astrophysical background processes is likely to be the 
right approach. Furthermore, any inconsistencies in the hadronic generators 
to describe the (anti)proton and (anti)nucleon spectra automatically result in a 
shift of the po parameter. For instance, an underproduction of antiprotons with 
a generator compared to available data results in a higher pq value when try¬ 
ing to describe measured (anti)deuteron data. Moreover, (anti)neutron spectra 
are very challenging to access in typical particle physics detectors. Therefore, 
the typical approach is to assume that the antiproton and antineutron pro¬ 
duction cross-sections are equal, and hence it is even more challenging to tune 
hadronic generators to (anti)neutrons than to antiprotons. Any possible pos¬ 
sible isospin asymmetry is an additional component of the po value. On the 
quantum-theoretical level, it is also important to realize that the formation 
probability in the per-event simulation approach is taken to be exactly 100% 
if the Pq condition is met. This is unlikely to be true when considering, e.g., 
that (anti)proton and (anti)neutron need to have aligned spins to form a stable 
(anti)deuteron with spin 1. Therefore, trying to separate the actual coalescence 
process from other conditions that have to be met will be an important challenge 
for the understanding of (anti)deuteron production. 

3.1.4- Alternatives to the coalescence model 

Few alternatives to the coalescence model have been suggested. Recently,the 
authors of |128j discussed a model in which the combination of an antinucleon 
pair with center-of-mass frame momentum difference k = \pp — Pn\ into an 
antideuteron, pfi —>■ dX, is a random event with a probability given by the ratio 
of the cross section of the corresponding scattering process as a function of k, 
to a normalization cross section, crpn^p[x{k)/<J q. The normalization tJo is a free 
parameter to be fixed through calibration against experimental data, analogous 
to Pq in the coalescence model. The main consequence of this approach is that 
it shifts the typical value of k involved in the antideuteron formation: where 
the coalescence model restricts the phase space to k < pq, with pq typically of 
the order of 100 — 200 MeV, this model prefers values of GeV where resonant 
production through the delta-resonance can occur. In turn, this affects the 
antideuteron spectrum, and the authors observe signihcantly improved fits to 
antideuteron production in recent p-p-data, from ALICE |129) . 
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3.2. Collider measurements of antideuteron production 

In order to make progress in the understanding of antideuteron formation 
and the prediction of the primary and secondary antideuteron fluxes, more ex¬ 
perimental data and a better understanding of the physics and modeling of 
antideuteron formation are needed. Recently, new experimental data that can 
be used for the determination of the coalescence momentum have become avail¬ 
able, and further determinations of the coalescence momentum will likely also 
be possible In the near future. 

The BaBar experiment measured the production rate of antideuterons in 
T(n«S') decays (with n = 1,2,3) (Sec. 3.2.1) as well as In continuum electron- 
positron annihilations to quarks at -^/s = 10.58 GeV, which can be used to tune 
Po in these processes. This determination, especially the one from e“'"e“ -A 
qq, Is particularly relevant for the calculation of the antideuteron flux from 
annihilations of light dark matter particles Into quark-antiquark pairs. 

The ALICE experiment at the LHC measured the production rate of an¬ 
tideuterons in p-p collisions at center-of-mass energies of ^/s = 7, 8TeV and 
will measure at 13TeV and 14TeV In the future (Sec. |3.2T2 ). Although at a 
different center-of-mass energy, this is precisely the main process participating 
in the antideuteron production by cosmic-ray spallations, hence the relevance of 
this process In calculating the secondary antideuteron flux. Collisions involving 
heavy ions, such as p-Pb and Pb-Pb collisions, are less relevant for studies of 
cosmic antideuterons, due to the different dynamics involved In these processes 
compared to dark matter annihilations, decays or cosmic-ray collisions. 

To reduce the systematic uncertainties in the calculation of the secondary 
antideuteron flux and to correctly Interpret a hypothetical future detection of 
cosmic antideuterons, an experiment measuring antideuteron production in p-p 
collisions at low center-of-mass energies, l.e. ^/s « 10 GeV, would be of ut¬ 
most Importance. This is precisely the main process and the most relevant en¬ 
ergy for the production of antideuterons by cosmic-ray spallations, constituting 
the background for dark matter searches. In addition, studies of antideuteron 
production in processes like p-C are important for constraining instrumental 
backgrounds and reducing systematic effects. The operational fixed target ex¬ 
periment NA61/SHINE is ideally suited for these tasks (Sec. 3.2.3). 

In addition, the upcoming PANDA experiment at the FAIR collider in Darm¬ 
stadt, Germany will probe the production of antideuterons in interactions of 
antiprotons with different fixed targets (Sec. 3.2.4). This is also an important 
process for the prediction of the astrophysical antideuteron background and will 
be discussed more In Sec. |4.1.1[ Using different beam energies and target mate¬ 
rials the cross section can be studied as a function of increasing energy starting 
slightly below the antideuteron production threshold. 

In the longer term, a high luminosity electron-positron linear collider with 
center-of-mass energies y/s « 10-500 GeV would be of great interest to investi¬ 
gate the antideuteron production in a single channel at different center-of-mass 
energies. Even reanalyses of LEP data could be helpful here. In particular, such 
study will allow to compare and tune different Monte Carlo generators and to 
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elucidate whether po is a constant or on the contrary scales with the center- 
of-mass energy. Moreover, the final state qq corresponds to the final state in 
the annihilation or decay in many well-motivated dark matter scenarios. Be¬ 
sides, the proposed super-B factory would deliver much more precise data on 
antideuteron production in T decays, thanks to the much larger luminosity 
compared to BaBar and BELLE. In particular, in combination with data on d 
production in the e“'"e“ continuum, it could help to test the universality of the 
coalescence momentum with respect to the underlying hard process. 

The following subsections will elaborate more on the different experimental 
aspects. 

3.2.1. BaBar 

Recently, the BaBar collaboration reported precise measurements of an¬ 
tideuteron production in hadronic decays of T resonances {T{nS), {n = 1,2, 3)), 
as well as the first measurement of antideuteron production in continuum (non¬ 
resonant) annihilation at ^/s ~ 10.58 GeV |123| . A brief description of the 
measurement is provided below as an illustration of the techniques employed to 
measure antideuteron production rates at e“'"e“ colliders. 

The measurements are performed using data collected with the BaBar de¬ 
tector at the asymmetric-energy PEP-II electron-positron collider at SLAG. 
The data used are collected running at e+e“ center-of-mass (CM) energies cor¬ 
responding to mx( 2 S )5 wx( 3 S) and mx( 4 S), along with data collected 40MeV 
below each resonance. The asymmetric energy of the collisions means that the 
center-of-mass (CM) frame is boosted in the laboratory (detector) frame by a 
factor P'y « 0.56. 

Deuteron and antideuteron candidates are selected from well-reconstructed 
charged particle candidates with lab-frame momentum between 0.5 and 1.5 GeV 
which are well-contained within the acceptance of the detector’s 40-layer mul¬ 
ticell drift chamber (DCH) (the BaBj^r detector is described in detail in |130j L 
Thanks to the boost of the collision frame, this corresponds to a wide range 
of deuteron CM momenta. This momentum range is below the Cherenkov 
threshold for antideuterons in the quartz radiator bars of the Cherenkov par¬ 
ticle identification system (DIRC) used at higher momenta, thus signals in the 
DIRC are used to veto lighter species. Finally, candidates are required to leave 
a well-reconstructed pattern of ionization in the DCH, with at least 24 good 
measurements of deposited charge sampled by the system. 

The resulting candidates are binned by momentum in the center-of-mass 
frame of the colliding e“'"e“ beams. The candidates are weighted by the inverse 
efficiency for triggering, detection and reconstruction, as well as by a factor 
which corrects for the fiducial acceptance of the detector at a given deuteron 
CM momentum. For each candidate, specific ionization {AE/dx) measurements 
performed in the DCH and in the 5-layer Silicon Vertex Tracker (SVT) are cor¬ 
rected for the differing materials and combined in a weighted average that also 
accounts for the expected resolution of the two independent measurements. The 
variable of interest is the distribution of normalized residuals of these averaged 
measurements with respect to the most probable value given the deuteron mass 
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Figure 7: Measured antideuteron differential spectra in e'*"e“ ^ g'g at a CM energy of ^/s ~ 
10.58 GeV. The points with inner (black) error bars give the measurements and their associated 
statistical uncertainties, the outer (red) error bars give the quadratic sum of the statistical 
and systematic uncertainties, and the dashed (blue) curve shows a fit to determine the total 
rate, from [I23]. 

hypothesis. For real deuterons and antideuterons, this distribution is nearly 
Gaussian and peaks at zero. The distribution for other species lies mostly at 
large negative values and falls rapidly near the signal region. 

A weighted unbinned simultaneous maximum-likelihood fit is used to ex¬ 
tract the yield of antideuterons in each bin. Matter deuterons are produced 
at a much higher rate via interactions of primary particles with detector and 
beampipe materials. These provide a pure high statistics sample which allows 
the parameters of the signal distribution to be determined in the simultane¬ 
ous fit by data rather than simulation. This fit is performed separately for 
BaBar's T{nS), (n = 2,3) datasets, a subset of the T{2S) dataset with well- 
reconstructed T(2S') — 7 . T(lS') 7 r+ 7 r“ candidates, and in the large ^/s = tox( 4 S) 
dataset. The latter is used to measure production in continuum, as direct T(4S') 
or B meson decay to antideuterons is expected to be negligibly small. This as¬ 
sumption is validated by cross-checks against data collected at a CM energy 
40MeV below mx(is)- 

As a result of the event-by-event weighting, the fits provide directly the dif¬ 
ferential cross section (or differential decay rate) for antideuteron production. 
For T(2,35') decay, below-resonance datasets are fitted simultaneously to sub¬ 
tract any continuum contribution. For T(IS'), the background subtraction is 
performed by using sideband regions in the variable 

^recoil — 's/ (Abeam (Pbeam Ptttt)^- (4) 

For the continuum measurement, the below-resonance yields are added to the 
yields in the on-resonance dataset. 

The results for continuum e^e“ annihilation are shown in Fig.[^ The results 
are in good agreement with previous measurements of antideuteron production 
in T(15') and T(2S') decay, and represent the first measurement of production 
in continuum e+e” annihilation with CT(e'''e“ —>■ dX) = (9.63 ± fb 

at ^/s « 10.58 GeV. Normalizing the continuum cross section to crg+e-_,,hadrons 5 
one finds a suppression relative to production in hadronic T decay of (3.01 ± 
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Figure 8: Deuteron spectrum for p-p collisions at y/s = 7TeV. The dotted line represents the 
fit with the Levy-Tsallis function. 


0.13lo;|j) • 10 as expected by the higher baryon production observed in gluon 
fragmentation relative to quark fragmentation. 

3.2.2. ALICE 

The results presented in this section are obtained from A Large Ion Collider 
Experiment (ALICE) at the LHC. Its performance and the description of its 
various subsystems are discussed in |13I] . Due to its unique particle identifi¬ 
cation capabilities ALICE is ideally suited to measure rarely produced nuclei. 
The collision energy reached at the LHC provides the opportunity to measure 
nuclei and the corresponding antiparticles in unprecedented abundances, al¬ 
though the measurement is challenging as the production probability decreases 
with increasing particle mass. The results presented here are obtained with 
the data of p-p collisions at -y/s = 7TeV recorded in 2011, Pb-Pb collisions at 
yfsNN = 2.76 TeV recorded in 2010 and p-Pb collisions at ^/sNN = 5.02 TeV 
recorded at the beginning of 2013. The following data analysis and discussion 
focuses on deuteron production. Symmetry arguments can be used to transfer 
this discussion to the production mechanism of antideuterons. 

Nuclei and antinuclei are identified over a wide transverse momentum {pt) 
range using the combined information of the specific energy loss (dE/da;) mea¬ 
surement in the Time Projection Chamber (TPC) |131j and the velocity mea¬ 
sured by the Time Of Flight detector (TOE) |132) . The measured energy loss 
signal of a track in the TPC is required to be within a three standard deviation 
region around the expected value for a given mass hypothesis described by the 
Bethe-Bloch formula |133j . Since these curves start to overlap at high momenta, 
the velocity measurement with the TOE is used in addition to allow for clear 
identification at pr > 0.8GeV/c. Fig. shows the efficiency and acceptance 
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Figure 9: Deuteron-to-proton ratio as a function of charged-particle multiplicity at midrapidity 
for different colliding systems at different energies. 


corrected deuteron spectra for p-p collisions at a/s = 7 TeV. The measurement 
of deuteron transverse momentum spectra has been also performed using the 
data collected during the p-Pb run. The deuteron spectra show a hardening 
with increasing multiplicity/centrality that is qualitatively similar to the pro¬ 
ton spectra. In heavy ion collisions this hardening is commonly attributed to 
radial flow |134j . 

Fig. [9] shows the deuteron-to-proton ratio for minimum bias p-p collisions 
and as a function of the multiplicity for p-Pb and Pb-Pb. The ratio rises with 
multiplicity until a saturation within errors in Pb-Pb collisions is reached. The 
transition between the different collision systems suggests that the d/p ratio 
is in part determined by the event multiplicity, at least for smaller systems 
{N < 100). 

The production spectra of light nuclei can be understood based on the coa¬ 
lescence approach assuming that deuterons (and other light nuclei) are produced 
by protons and neutrons that are close in phase space. In the most naive picture, 
this would lead to an increased deuteron production for higher nucleon multi¬ 
plicities. The increase of the deuteron-to-proton ratio with the charged particle 
multiplicity in Fig. is consistent with this picture for small colliding systems 
like p-p and p-Pb. However, the d/p ratio for Pb-Pb collisions is constant with 
increasing centrality, although the nucleon multiplicity increases. 

A possible explanation of the Pb-Pb results can be that the increasing nu¬ 
cleon multiplicity is counterbalanced by the increasing source volume, leading to 
a constant nucleon density. This is also consistent with the rising d/p ratio with 
multiplicity in p-Pb collisions, if the effect of the increasing nucleon multiplicity 
dominates over the effect of the increasing source volume. 

In addition to the coalescence model presented in Sec. |3.1| the production 
mechanism of antinuclei in heavy ion collision can also be discussed in a thermal 
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model |135L 1136] . The production yields depend exponentially on the chemical 
freeze-out temperature Tchem and the mass m: 



(5) 


Due to their large masses the abundance of nuclei is very sensitive to Tchem- An 
important question is whether the nuclei are produced at the chemical freeze- 
out or at a later stage via coalescence. Besides the constant d/p ratio a key 
observation is that in Pb-Pb collisions the absolute production yields {dN/dy) 
of light nuclei are in good agreement with thermal model calculations, as shown 
in [137] . On the other hand the highest d/p ratio obtained in p-p collisions is 
about half the value predicted by a thermal model in central AA collisions, dis¬ 
favoring the statistical description based on thermal equilibrium in p-p collisions 
at the LHC energies. Further studies are needed to establish whether the fast 
expansion conserves the particle ratios and which additional conditions in the 
coalescence model are required to describe the constant particle ratio in Pb-Pb. 

3.2.3. NA61/SHINE 

The NA61/SHINE (SPS Heavy Ion and Neutrino Experiment) experimental 
facility at the Super Proton Synchrotron (SPS) at CERN is used to study col¬ 
lisions of protons and heavy nuclei with fixed targets in a wide incident beam 
momentum range (e.g., 13-158 GeV/c protons on a hydrogen target) |138j . An- 
tideuterons are also produced in these interactions and NA61/SHINE is cur¬ 
rently the best operational experiment to study antideuteron production from 
the threshold on. Different subdetectors are needed for the particle identifica¬ 
tion. Numerous of these components were inherited from its predecessor, the 
NA49 experiment [139] . The main upgrades from NA49 to NA61/SHINE in¬ 
clude a forward time-of-flight detector wall, improved time projection chamber 
readout, and five times better energy resolution with the projectile spectator 
detector. 

A set of scintillation and Cherenkov counters as well as beam position de¬ 
tectors upstream of the spectrometer provide timing reference, identification 
and position measurements of incoming beam particles. A trigger scintillator 
counter placed downstream of the target is used to select events with collisions 
in the target area. The main tracking devices of the spectrometer are large 
volume time projection chambers (TPC). Two of them are located in the mag¬ 
netic fields of two super-conducting dipole magnets with a maximum combined 
bending power of 9 Tm, which corresponds to about 1.5 T in the first and 1.1 T 
in the second magnet. Two large TPCs are positioned downstream of the mag¬ 
nets symmetrically to the beam line. The fifth small TPC is placed between 
the two vertex time projection chambers on the beam line. It closes the gap be¬ 
tween the beam axis and the sensitive volumes of the other TPCs. The particle 
identification capability of the TPCs is based on measurements of the specific 
energy loss, dE/dx, and is augmented by velocity measurements using time- 
of-flight detectors {at = llOps). The high resolution forward calorimeter, the 
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projectile spectator detector, measures energy flow around the beam direction, 
which in nucleus-nucleus collisions is primarily given by the projectile spectators. 
NA61/SHINE uses various solid nuclear targets and a liquid hydrogen target. 
The targets are positioned about 80 cm upstream of the sensitive volume of the 
first vertex time projection chamber. 

NA61/SHINE can be used to improve the parameters of the underlying 
(anti)deuteron coalescence model. Proton-proton interactions with incident mo¬ 
mentum between 13 and 158GeV/c were already recorded in 2009 and 2011 
and data with higher energies will be taken in 2015. The low-energy cosmic-ray 
antideuteron background flux (< IGeV/n) is dominated by products of pro¬ 
ton interactions with interstellar medium hydrogen at somewhat higher kinetic 
energies (2-5 GeV/n), which is then shifted towards lower energies due to prop¬ 
agation through the Galaxy [^. Galculations predict that the most relevant 
cosmic-ray production happens between 40 and 400 GeV/c. Therefore, the SPS 
energies from 9 to 400 GeV/c are ideally suited for a study relevant to the an¬ 
tideuteron search [4^. In addition, NA61/SHINE scans different ion collisions 
and different energies like p-G, tt-G, Be-Be, Ar-Sc. Especially the production of 
antideuterons in p-G interactions is important to understand systematic exper¬ 
imental effects due to antideuteron production inside experiments. 

The NA61/SHINE (anti)deuteron production can be studied as a function 
of rapidity and transverse momentum and used to calibrate the coalescence 
model. NA61/SHINE can also be used to test the reasonable assumption that 
the merging of pn and ph pairs should proceed similarly for nucleons and antin¬ 
ucleons. For this purpose, it will be important to also measure deuterons 
with NA61/SHINE at the same time. However, as discussed at the end of 
Sec. |3.1.3[ only specific deuteron channels allow to draw direct conclusions 
for antideuterons. For instance, the antideuteron production cross section in 
pp —> dnp should be the same as for deuteron production in pp —> dnp. The 
advantage of measuring pp —>■ dnp in comparison to pp —> dX is that the 
cross-section is larger, because only four instead of six (anti)nucleons have to 
be produced. Therefore, using the pp —>■ dnp channel can already answer the 
question of how (anti)deuteron coalescence depends on the available energy. 
Another important cross-check for the Monte Garlo generators will be the mea¬ 
surement of the yield of antiprotons with the same data. The Pb-Pb data of the 
predecessor of NA61/SHINE, NA49, were already analyzed for antideuterons 
and successfully demonstrated that antideuterons can be identified in the ex¬ 
perimental setup |140) . However, it is theoretically not clear how to apply the 
event-by-event coalescence model to proton-nucleus or nucleus-nucleus collisions 
because the physics of the creation of secondary nucleons is very different in this 
case. 

3.2.4- The PANDA fixed target experiment 

The upcoming PANDA experiment at the FAIR collider in Darmstadt, Ger¬ 
many will study the interactions of a 1.5-15 GeV/c antiproton beam with dif¬ 
ferent fixed targets (e.g., hydrogen, deuterium, or even gold), and thus will 
also probe antideuteron physics [T 4 n[T 42 ] . It will reach luminosities of up to 
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2 • corresponding to a pp-annihilation rate of 2 • 10^ s“^. The 

PANDA detector is a multi-purpose particle physics detector with 47r accep¬ 
tance. The detector is divided into two main components, namely the target 
and the forward spectrometer. 

The target spectrometer surrounds the interaction point and is split by the 
target injection pipe into two half shells. The typical target will be gaseous 
hydrogen or other heavier gases as well as solid targets for hypernuclei studies. 
The detector concept follows a layer structure. The interaction point is sur¬ 
rounded by a microvertex detector including radiation hard silicon pixels. The 
next layer is the central tracker and will consist of a 24-layer straw tube tracker. 
The forward direction will be covered by three layers of GEM tracking stations. 
Particle identification of slow particles (< IGeV/c) for the barrel region will 
be done by a time of flight detector with a timing resolution of 50-100 ps and 
the specihc energy loss in the central tracker. Faster particles are identified 
by Cherenkov detectors for the barrel and forward endcap. The next layer is 
formed by the electromagnetic calorimeter made out of lead tungstate with a 
short radiation length and Moliere radius to assure a compact design. All the 
former mentioned subdetectors are contained in a 2 T magnetic held provided 
by a superconducting solenoid coil. Outside of the magnetic held as part of the 
magnetic return yoke a hne segmentation with aluminum drift tubes serves as 
a range tracking system for muons and pions. In the forward direction particles 
are dehected by a dipole magnet. The particle identihcation of the forward spec¬ 
trometer is similar to the target spectrometer. The forward tracker measures 
the dehection in the magnetic held with three pairs of tracking detectors based 
on straw tubes. Particle identihcation will be done by a ring image Cherenkov 
detector and a time of hight system followed by a lead-scintillator sandwich 
electromagnetic calorimeter. Forward muons are again detected with a range 
tracking system laid out for higher momenta. At the very end of the line, a 
luminosity detector measures elastically scattered antiprotons with four silicon 
detector layers. 

Using different beam energies and target materials the cross section can 
be studied as a function of increasing energy starting slightly below the an- 
tideuteron production threshold with hydrogen. Furthermore, by analyzing 
antiproton-deuterium interactions the inelastic scattering cross section of an- 
tideuterons with matter can be studied assuming CPT-symmetry. 

3.3. The coalescence model for antihelium 

The yield of antinuclei produced by dark matter annihilation or decay pro¬ 
cesses decreases rapidly with the antinucleus mass number. As a rough estimate, 
one can assume that the probability to have coalescence is reduced by a factor 
0(10“^) for each additional antinucleon that participates in the process. This 
is why only antihelium formed by three antinucleons (^He), will be taken into 
account. The production rate of ^He is even more suppressed and will not be dis¬ 
cussed further. Artihcially-produced antihelium-3 was discovered in 1971 |143) 
and antihelium-4 in 2011 |144j . 
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The coalescence process that leads to the production of ^He can involve either 
two antiprotons and one antineutron (in which case the antihelium is produced 
directly) or two antineutrons and one antiproton (in which case the antihelium 
comes from the decay of the antitritium, a phenomenon that occurs on a time 
scale that is much shorter than those characterizing Galactic propagation and 
can therefore be considered as instantaneous for our purposes). Coalescence 
in the (p, p, n) channel is expected to be suppressed by Coulomb repulsion due 
to the presence of two antiprotons [55]. If this suppression is not taken into 
account, the two different coalescence channels instead have practically the same 
yield mum]. 

Two Monte Carlo coalescence models for antihelium have been proposed. 
The first requires that the relative momenta of all sets of the three nucleons 
are smaller than the coalescence momentum po m, while the second requires 
that each of the relative momenta lie within a “minimal bounding momentum 
sphere” with diameter po m- The flux predicted using each model can vary 
by, at most, 15% [112j . 

Laboratory measurements of the antihelium production rate are extremely 
scarce and related to processes that are quite different from dark matter annihi¬ 
lation, such as proton-nucleus [1451114611147] or nucleus-nucleus collisions [148] . 
Therefore, the value of the coalescence momentum used for dark matter signal 
predictions is highly uncertain. Ref. [Ill] adopted the same fiducial coalescence 
momentum determined for the antideuteron case, i.e. 195 MeV; Ref. m in¬ 
stead rescaled the coalescence momentum derived for the antideuteron case by 
a suitable factor. This factor can either be determined by assuming that the 
coalescence momentum is proportional to the square root of the binding energy 
B of the bound state: 


3 

Po 


He 



(0.357 ± 0.059) GeV 


( 6 ) 


or by assuming a proportionality between Pq^® and Pq inferred by comparing 
theoretical expectations with the experimental measurements of |148j : 


Po“® = 1.28p(( = (0.246 ± 0.038) GeV. (7) 

To ensure that the p and n originate from the same production vertex, as 
has been discussed for the antideuteron case, all decays of long-lived particles 
in the Monte Carlo event generator are turned off. 


4. Antideuteron propagation 

After coalescence, antideuterons must propagate through the Galactic mag¬ 
netic fields, plasma currents, and the interstellar medium to reach our solar 
system. They can then be deflected by the solar magnetic field or suffer adi¬ 
abatic energy losses in the solar wind. Finally, antideuterons can be deflected 
away from balloon-borne or satellite detectors by the Earth’s geomagnetic field 
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and/or interact with Earth’s atmosphere. Each of these processes impacts the 
predicted antideuteron flux seen by experiments, typically calculated at the top 
of atmosphere (TOA). Galactic propagation is an important uncertainty in pre¬ 
dicting antideuteron fluxes at Earth, with well-motivated choices of transport 
model parameters predicting antideuteron fluxes that can differ by an order of 
magnitude. Recent positron data exclude the MIN Galactic propagation model, 
which predicts the lowest antideuteron flux levels at Earth, supporting higher 
antideuteron flux predictions that could be detected by GAPS or AMS-02. 

The propagation of antideuterons in the Galactic environment is detailed in 
Sec. |4.1| Transport in the solar environment, geomagnetic field, and atmosphere 
is discussed in Sec. |4.2p^ and followed by a brief introduction of antideuteron 
interactions with detectors (Sec. |4.5[ ). The section concludes by elaborating on 
a possible extragalactic origin of antideuterons (Sec. 4.6). 

/.I. Transport in the Galactic environment 
4-1.1. Cosmic-ray transport equation 

Once produced, antinuclei propagate inside the halo of the Milky Way, where 
they can be deflected by the Galactic magnetic fields and local plasma currents, 
lose energy through interactions with the interstellar medium (ISM) and mag¬ 
netic fields, and be destroyed by fission or annihilation. Although the magnetic 
fields are significantly turbulent, cosmic-ray transport can still be modeled by 
spatial diffusion, which plays a major role above about 10 GeV, and convec¬ 
tion, which contributes primarily at lower energies. Gosmic-ray fluxes are thus 
determined by a transport equation as given, e.g., in |149j . 



The left-hand side describes spatial diffusion [K) and convection (V/), and 
the first and second order energy transport terms (/o and So)- The right-hand 
side corresponds to the source term {qsrc{T,E)), and the sink term describing 
destruction in the interstellar medium (ISM) (Tdst)- 

The diffusion coefficient K depends a priori on the location r and energy E 
of the particles. It is related to the power spectrum of the magnetic inhomo¬ 
geneities, which is poorly known. Although many forms have been proposed, it 
is often assumed to depend only on energy as [Ml EH]: 



(9) 


where the normalization Kq and the exponent S are assumed to be constant 
parameters, and (3 and TZ = pc/{Ze) denote the cosmic-ray velocity and rigidity, 
respectively. 

At energies below about 10 GeV, additional processes come into play. Antin¬ 
uclei are swept by the convection of the local plasma, which results from stellar 
winds, and drift with a velocity V/ upwards and downwards with respect to 
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the Galactic Disk. The conservation of the cosmic-ray current in energy space 
yields the first and second order energy transport terms. The first order term, 
fo{r,E), corresponds to the sum of four processes: ionization, Coulomb, and 
adiabatic losses, and first-order reacceleration. Ionization losses take place in 
the neutral ISM, while Coulomb losses are dominated by scattering off thermal 
electrons [15211153| in the completely ionized plasma. The spatial dependence 
of these two terms is thus encoded in the distribution of the neutral and ionized 
gas. Adiabatic losses are due to the expanding stellar winds, and their spatial 
dependence is thus related to the gradient of 14- 

Both the contribution to /o(r, A) from reacceleration and the second-order 
term So(d E) originate from the motion of knots in the turbulent Galactic mag¬ 
netic fields. In addition to being responsible for spatial diffusion, these Alfvenic 
waves also cause energy drift and reacceleration. A minimal reacceleration 
scheme is well-motivated m and allows us to calculate the /o and So coef¬ 
ficients. Similar, albeit more empirical, forms have also been used [15511156] . In 
all these models, the strength of the reacceleration is mediated via the Alfvenic 
speed Va of the scatterers. 

On the right-hand side of Eq. ([^, qsrc{r, E) denotes the rate at which antin¬ 
uclei are produced. Primary antinuclei refer to those produced by the annihi¬ 
lation of dark matter species in the Milky Way halo. For antideuterons, this 
production rate can be written as [H] US] : 

Ej) = - {av) \ dark matter annihilation, 

2 dAj V moM / 

qy^{f,Eg) = dark matter decay. 

tdm dEj todm 

Here, {av) is the thermally-averaged dark matter pair-annihilation cross sec¬ 
tion, Tdm is the dark matter lifetime, dN^/dE^ is the injection spectrum of an¬ 
tideuterons produced by individual annihilations or decays, pdm( 4) is the dark 
matter density at location T, and todm is the dark matter particle mass. Sec¬ 
ondary antinuclei refer to those produced via interactions of high-energy cosmic- 
ray protons and helium nuclei on the ISM. These are the irreducible background 
for any dark matter annihilation signature. In the case of antideuterons pro¬ 
duced by cosmic-ray protons impinging on hydrogen atoms at rest, the secondary 
production rate can be written as nails]: 


( 10 ) 

( 11 ) 


|•+oo 

qf^{r,Es)= -Jl^0^{Ep^Es)n^{r)vpNj,{r,Ej,)dEj,, ( 12 ) 

Jeo dAj 

where uh is the hydrogen density in the ISM, Vp is the velocity of the protons, 
and Np is the differential proton density, Np{f,E) = dup/dEp. The differen¬ 
tial production cross section, dap^^^x computed using the coalescence 
scheme described in Sec. and can be adapted to account for production via 
helium in the ISM and in the cosmic radiation field. 
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Also cosmic-ray antiproton collisions with the ISM can produce antineu¬ 
trons, which can then merge with cosmic-ray antiprotons into antideuterons. 
Production from proton and antiproton collisions with the interstellar medium 
have different kinematics. For the proton production enough energy for at least 
six nucleons must be present while for antiprotons only four are needed because 
the incident antiproton can coalesce with an antineutron collision product. For 
the case of impinging protons and antiprotons on interstellar medium hydro¬ 
gen at rest, the corresponding kinetic energy production thresholds are 17 GeV 
and 7 GeV, respectively. Therefore, the two dominant production processes are 
pp —>■ dX and, although the antiproton flux is strongly suppressed compared 
to protons, pp dX. The antiprotonic production is especially for the lower 
energy part of the spectrum non-negligible |49j . 

As antinuclei propagate inside the Galactic Disk, they can be destroyed 
through fission or annihilation with the ISM. In the case of antideuterons, the 
annihilation rate F^st is given by m- 

^dst (riH -f 4^/^nHe) '^d ^ine {dp ^ V), (13) 

where the hydrogen and helium densities in the ISM are assumed homo¬ 
geneous with TiH = 0.9cm“^ and riHe = 0.1 cm“^, uj is the velocity of the 
antideuteron, and the cross section cTine accounts for inelastic interactions of 
antideuterons on the ISM. To calculate this cross section, which has not yet 
been measured, one can follow the assumptions outlined in |56j . See also the 
discussion in [6^ . 

Antideuterons rarely survive inelastic collisions. However, they can some¬ 
times undergo inelastic, but non-annihilating collisions, and transfer enough 
energy to excite the target proton as a A resonance. As they typically lose 
a substantial fraction of their initial energy, these antideuterons are effectively 
re-injected into the total flux with a degraded momentum. This mechanism 
redistributes antideuterons towards lower energies and flattens their spectrum, 
as first remarked by [13 II57j . The antideuterons that result from these inelastic 
collisions are referred to as tertiary antideuterons, with a production rate given 
by m - 


/ + 00 Aq-- 

nH(rl u'-iV(f, A' 




^dp^dx {Ed) riR {r) N{r, Eg). 


(14) 


In this expression, the differential cross section for inelastic and non-annihilating 
interactions is often approximated by [157] : 


(15) 


do-dp^dx ^ ^ ^dp^dx {Eg) 

-^^(Ed^Eg) = y, > 

where Tg denotes the initial antideuteron kinetic energy. To account for the 
tertiary production rate Eq. (141 in the transport equation Eq. ([^, one needs to 


31 










subtract the cross section for inelastic, non-annihilating interactions from (Tine 
in the calculation of the annihilation rate Fdst- 

Once the density N is derived from Eq. ([^, the interstellar flux at the solar 
system is given by |42) : 

^d{E) = ^ n^{r = rQ,z = 0,E) (16) 

_ ^ ^ - (av) —- for dark matter annihilation, 

^TrymuMj ^ 

(17) 

_ h Pq i^d®'=(£') for dark matter decay, (18) 

47r todm tdm aEg 

where the propagation function R^{E) encodes all the dependence on astro- 
physical parameters and is determined numerically. The specihc form of R^{E) 
can be found in |42) . 

4-.1.2. Two-zone diffusion model 

A full numerical treatment is generally required to solve the transport equa¬ 
tion Eq. ([^, as described, e.g., in |158) . However, analytical (or semi-analytical) 
solutions may be derived assuming a simplified description of the spatial depen¬ 
dence of some parameters. The two-zone diffusion model [15911160| . based on 
the description of the Galaxy as a thin gaseous disk embedded in a thick diffu¬ 
sive halo, has proven to be successful in reproducing the nuclear [1611115611162] , 
antiproton [163] . and radioactive isotopes [162] data. It also allows us to treat 
contributions from dark matter (or other exotic) sources located in the magnetic 
halo [16411165L1166] . This model has been extensively detailed in [15611167L1168] . 

In the two-zone model, the Galaxy is defined as a cylinder with a magnetic 
halo of half-height z = L and radius r = R. The ISM and the nuclei accelerators 
are contained in a thin-disk of half-height h L. The halo thickness L is a 
free parameter of the model, and h and R are set to 100 pc and 20 kpc, respec¬ 
tively. The diffusion coefficient Eq. is assumed to be the same throughout 
the Galactic magnetic halo, and the convective wind is assumed to be of con¬ 
stant magnitude directed outwards perpendicular to the Galactic Plane with 
Vc = VcBz. The reacceleration strength, mediated by the Alfven velocity 14, is 
confined to the thin disk. The first and second order terms /□ and Sq in Eq. (§ 
follow the formulation given in [156] . 

In the framework of the two-zone diffusion model, the details of Galactic 
propagation depend entirely on the values that are assigned to the parameters 
(L, 14, Eq, 6, 14). These parameters are usually constrained by measurements 
of cosmic-ray observables such as the primary fluxes and secondary-to-primary 
ratios (in particular B/G). However, fitting to existing B/G data does not lift the 
strong degeneracy between these hve transport parameters. In particular, many 
sets of these parameters lead to the same B/G ratio and the same secondary 
(standard) antiproton flux [163] . This degeneracy is broken for sources located 
in the diffusive halo, leading to large astrophysical uncertainties for the relevant 
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fluxes |164LI169LI170| . In particular, a large uncertainty for dark matter searches 
stems from the degeneracy between the normalization of the diffusion coefficient 
Kq and the halo size L. Indeed, stable secondary-to-primary ratios mostly 
constrain the cosmic-ray escape time oc L/Kq. However, L has a strong impact 
on dark matter signal predictions since it is proportional to the amount of cosmic 
rays induced by dark matter annihilation or decay within the diffusive Galactic 
Halo that contribute to the local cosmic-ray flux. 

As benchmark scenarios, the three sets of parameters labeled as MIN, MED 
and MAX defined in |164) are adopted in this review. A large range of L values 
can be found in the literature. The values L = 1 and 15 kpc are used in the 
MIN and MAX models, respectively. These values were proposed by |164j to 
bracket the theoretical uncertainties on dark matter signal predictions, relying 
on the boron-to-carbon (B/C) analysis performed in |156j . Fig. 10 shows the 
ratios of antideuteron fluxes from dark matter annihilation (todm = 100 GeV, 
bb channel) in the MAX and MIN propagation models with respect to the MED 
propagation model. For low energies the relative difference between MIN and 
MAX covers about one order of magnitude. Some correlation between the un¬ 
certainties for the propagation of astrophysical antideuterons and dark-matter 
induced antideuterons exists. However, this is not a dominant effect because 
the main uncertainty for the dark matter signal comes from the fact that it is 
produced in the halo and not, like the astrophysical background, in the Galactic 
disk. While the astrophysical background is nearly unchanged by the choice of 
halo size, a larger halo produces more dark-matter induced antideuterons. 


4-. 1.3. Recent experimental constraints on diffu.sion model parameters 

The Kq/L degeneracy can be broken by an observable that is sensitive to 
only one of the two parameters. In general, radioactive secondary-to-stable 
secondary ratios, such as ^°Be/®Be, have been used. The radioactive secondaries 
decay before they can reach the edge of the Galaxy and escape, and thus their 
modeling is independent of the Galactic diffusive halo size L. However, several 
studies have shown that this method is very sensitive to the modeling of the 
local ISM and strongly affected by the presence of a local under-density, known 
as a local bubble |171j . In addition, the lack of precise measurements over a 
sufficiently large energy range does not permit a precise estimation of the halo 
size. 

Indirect constraints on L can be obtained from calculations of diffuse Galac¬ 
tic y-rays m or radio |173L 1174) emissions, disfavoring very small (large) halo 
values L < 2 kpc (> 10 kpc). Nevertheless, predictions of these observables rely 
on additional model parameters with non-negligible systematics and correla¬ 
tions, for example line-of-sight integrals depending on the astrophysical source, 
interstellar medium, and/or magnetic field distributions. 

Recently, the authors of |175) demonstrated that low-energy secondary posi¬ 
trons can directly constrain two-zone diffusion models with small haloes and 
large diffusion indices 5, effectively excluding the MIN model. Secondary posi¬ 
trons are produced by inelastic scattering of primary cosmic rays off the ISM. 
Due to energy losses, the mean free path of secondary positrons is much smaller 
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Figure 10: Ratios of antideuteron fluxes from dark matter annihilation (moM = 100 GeV, bb 
channel) in the MAX and MIN propagation models with respect to the MED propagation 
model, based on 02]. 


than that of nuclei, decreasing the dependence on L. For small L models, bound¬ 
ary effects become important and the secondary positron flux scales roughly 
with 1/L. This excludes small L models which predict a higher than observed 
positron flux, and excludes high spectral indices through the spectral shape. 
Besides the availability of very precise positron data over a large energy range, 
this method is less sensitive to the modeling of the local interstellar medium 
compared to the standard approach. 

The exclusion curves shown for different spectral indices S (red solid lines) in 
Fig. are obtained by comparing 500 000 predicted secondary positron fluxes 
consistent with B/C constraints to the PAMELA data [176) . Models leading 
to fluxes in excess by more than three standard deviations with respect to the 
data are excluded. In the left panel, the nominal value of 520 MV for the solar 
modulation (Sec. 4.2) was adopted. In the right panel, conservative results are 
shown by (i) considering only the data above 2.38 GeV and (ii) increasing the 
solar modulation potential to 700 MV. In the same plots the B/C constraints of 
urn for the different slices of S are represented by green bands. The obtained 
bounds are almost orthogonal to the B/C constraints, and thereby uncorrelated. 
In the nominal and conservative case, the MIN propagation model, and more 
generally large diffusion indices (i5 > 0.8), are excluded. This is good news 
for cosmic-ray antideuteron searches, as the MED and MAX diffusion models 
predict at least an order of magnitude higher low-energy antideuteron fluxes, as 
seen in Eig. [T^ 
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Figure 11: Constraints on propagation parameters in the \og{KQ/L • 1 Myr/1 kpc)-L [kpc] 
plane. Lines are constraints from the positron flux (downward regions are excluded). Green 
filled contours are allowed by B/C data. The MIN and MED models of ESI] are indicated by 
a filled black circle. Left: Positron contours (from PAMELA data) for a realistic modulation 
level of ^ = 520 MV. Dashed lines correspond to the limits if the secondary positron prediction 
is increased by 30% to mimic a primary component. Right (super conservative): Same, but 
with a solar modulation of 700 MV and considering only data points above 2 GeV. Figures are 
taken from EZSl- 


Effect of dark matter density profiles 

Due to propagation effects, charged antideuterons in the relevant energy 
range cannot be traced back to the origin of dark matter annihilation or decay. 
Therefore, a proper determination of the effect induced by changing the dark 
matter profile cannot be separated from the size and properties of the diffusion 
zone. Nevertheless, a qualitative understanding can be obtained by realizing 
that the dark matter signals depend on the quantities: 


A = p{r)^dV for dark matter annihilation, (19) 

Jv 

D= p{r)dV for dark matter decay, (20) 

Jv 

with V being a cylinder with radius and height from the propagation model 
under study and p(r) the dark matter density. D is quite insensitive to the 
halo profile, while A instead has a more pronounced dependence. The reason is 
that for a cuspy profile most of the annihilations occur inside the diffusion zone, 
whereas for a cored profile, one can also have annihilations outside. In the latter 
case, some of the antimatter particles produced in the annihilation can escape 
into intergalactic space, and hence cause some differences in the fluxes at the 
Earth. This difference is not very significant, since the dark matter distribution 
is practically the same everywhere inside the diffusion zone except for the small 
region around the Galactic Center. Ref. [164] explicitly investigated the effect 
of dark matter profile choice together with confinement zone and transport 
parameters for dark matter annihilating into antiprotons and found a typical 
difference of about 20% or less for small confinement regions and a maximum 
difference of 70% for very large confinement volumes. For dark matter decay the 
uncertainty from the choice of the halo profile is at the level of 10% or less |177] . 

In addition to the more involved combination of dark matter profile and 
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propagation, the local dark matter density po just acts as a dark matter flux 
normalization parameter. A recent calculation urn finds a local dark matter 
density value of po = 0-42 with an error at the 10%-level, which results in a 
dark matter profile uncertainty for dark matter decay of 10% and dark matter 
annihilation of 20%. 


4 . 1 . 5 . Galactic transport of antihelium 

Galactic transport of antihelium nuclei can be modeled similarly to the an- 
tideuteron case, with suitable modifications concerning interactions with the 
ISM. The most relevant difference is in the cross sections for annihilating colli¬ 
sions between the antihelium nucleus and the hydrogen and helium nuclei of the 
ISM. These cross sections can be modeled with the parameterizations provided 
in Table 4.5 of |179j . Similar to the antideuteron flux, the antihelium interstellar 
flux can be expressed in terms of the propagation function R as in Eq. (181. A 
parameterization of this function for several choices of the dark matter density 
profile is provided in m- 


4 . 2 . Transport in the solar environment 

An important effect for cosmic-ray measurements in the GeV energy range 
is the modulation by the magnetic held of the Sun that depends on the solar 
cycle. When the antinuclei reach the heliosphere, they experience diffusion due 
to inhomogeneities in the solar magnetic held (SMF), begin drifting along the 
held lines of the SMF, and suher adiabatic energy loss in the solar wind. Sev¬ 
eral approaches, including the simple force-held approximation and numerical 
simulations, have been discussed in the literature during the last decades. 


4 . 2 . 1 . Transport equation 

Similar to the galactic environment, the propagation in the heliospheric en¬ 
vironment is described by a transport equation [I80j : 

= -(Fsw + fTd) • V/ + V • (K • V/) + I (V • F,w) ^ , (21) 

where / is the antideuteron phase space density averaged over momentum 
directions, p is the antideuteron momentum, K is the (symmetrized) dihusion 
tensor, and the vectors 14w and Vd are the velocity of the solar wind and the 
velocity associated to antideuteron drifts, respectively. 


4 . 2 . 2 . Force-field approximation 

A simple solution to the propagation equation (21) is the force-held model 
approximation: 


Fn,od(A = Eus-\Z\e^) = F{r,Eus-\Z\e^) = F{oo,Eus)- S > (22) 

-^LIS 

where F(oo,Flis) is the hux in the local interstellar medium, Alis is the to¬ 
tal energy of the particle with mass toq and charge Ze and the effective solar 


36 











Figure 12: Effect of solar modulation for antideuterons as a function of solar modulated kinetic 
energy (notation as in Eg. |22[|. 


modulation parameter $ for all particle species, which is tuned to experimental 
data. The relative effect of solar modulation using the force-field approxima¬ 
tion for the antideuteron flux is shown in Fig. [^for three different choices of 
<f>, representing different periods of solar activity. For larger values of <f> the 
flux measured in the solar system experiences a stronger modulation. The solar 
modulation potential $ correlates with the sun activity, and thus low-energy 
antideuteron measurements are ideally carried out during solar minima. How¬ 
ever, the effect of varying solar modulation parameters within a typical range 
(from 300 to 700 MV at a kinetic energy of 100 MeV/n) is only at the 25% level. 
It is important to note that the parameter $ is known with about 10% accuracy 
from the measurements of other cosmic-ray species like protons |181j . Starting 
from about 20GeV/n, solar modulation becomes only a small effect. 


4-2.3. Numerical simulation 

An alternative approach to solving Eq. (21), 
in the HelioProp numerical code 


such as recently implemented 
consists in adopting a Monte Carlo 
treatment of heliospheric propagation, in which a sample of particles is injected 
at the modulation boundaries and traced to the position of the Earth. Here, 
the SMF is modeled as a large rotating spiral. Its geometry is characterized 
by the heliospheric current sheet, which is a surface separating magnetic field 
lines according to their polarity, with an exact shape depending on a parameter 
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Figure 13: Simulated number of days as a function of rigidity antiprotons are able to reach 
AMS-02 and GAPS through the geomagnetic field (horizontal lines indicate the means for 
different rigidity regions). 


known as the tilt angle. The mean-free-path of antideuterons in the direction 
parallel to the SMF can be expressed as [42] : 


Am = An 




1 GV 


B 


B, 


© 


-1 


(23) 


where B is the intensity of the SMF (i ?0 = 5 nT is the value of the mag¬ 
netic field at the Earth position [18311184] ) and the normalization Aq as well 
as the exponent 7 are assumed to be constant. In contrast to the force-field 
approximation, this model includes the important effect of Vd- As described in 
| 12 ]: the extent to which the solar modulation affects the antideuteron fluxes 
depends on four parameters: the polarity of the SMF, the tilt angle a, Aq, and 
7 . Ref. [12] showed that the difference between the more detailed approach and 
the force field approximation is small (« 20 %) compared to the uncertainties in 
the antideuteron coalescence model and Galactic propagation. However, a big 
advantage of the Monte Carlo approach is its ability to predict event-by-event 
fluctuations and to take into account spectral features of dark matter annihi¬ 
lation and background flux. This allows to use the solar activity level in the 
modeling when evaluating a particular dark matter model. 
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Figure 14: Attenuation of antiprotons in Earth’s atmosphere at 37 km altitude for an Antarctic 
balloon flight trajectory. 


4-3. Geomagnetic deflection 

The geomagnetic field plays a crucial role in the detection of charged cosmic 
rays. The geomagnetic field is roughly described by a tilted dipole field, which 
provides the strongest charged particle shielding at the equator and the weakest 
at the poles. Utilizing the Geant4-based PLANETOCOSMICS |185LI186] framework 
for backtracing charged particles through the Earth’s main magnetic field using 
the International Geomagnetic Reference Field |187j and the Tsyganenko 2004 
for the external magnetosphere |188) . Fig. 13 compares the number of days 


that antiprotons can reach AMS-02 (over five years of orbit) and GAPS (over 
three 35-day Antarctic flights) as a function of particle rigidity. About 10% 
of the measurement time AMS-02 is exposed to antideuterons of cosmic origin 
in the low-rigidity region (Sec. 5.2), which is most interesting for dark matter 
searches. The geomagnetic deflection is much lower above Antarctica and GAPS 
is sensitive 100% of the time. For AMS-02 the relative exposures increase in the 
higher rigidity regions to about 45% and 85%, respectively. These geomagnetic 
comparisons do not account for detector acceptances. It has also to be noted 
that solar magnetic disturbances influence the magnetic field in the vicinity 
of the Earth and can lead to considerable changes of the geomagnetic cutoff 
value [Ml nig, which have to be taken into account for the AMS-02 data 
analysis. 
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Atmospheric influence 


One more important effect is the influence of the atmosphere. This is di¬ 
rectly apparent when comparing the average amount of matter traversed by 
cosmic rays in the interstellar medium (6-10 g/cm^) to the overburden atmo¬ 
sphere (wGg/cm^). 


Fig. 


14 


shows results of simulations that have been car¬ 


ried out in the same Geant4-based approach as the geomagnetic simulations 
discussed above nmnng. These studies have two components: (i) the at¬ 
tenuation of cosmic-ray particles in the atmosphere and (ii) the atmospheric 
background production. Fig. compares the top-of-the-atmosphere kinetic 
energy of antiprotons to the average kinetic energy after reaching a typical bal¬ 
loon altitude of 37 km above Antarctica. The energy loss in the atmosphere is 
at the 10%-level, and thus the measured kinetic energy at balloon altitudes can 
be corrected to top-of-the-atmosphere kinetic energies without introducing large 
systematic uncertainties. Antideuterons at the same kinetic energy per nucleon 
lose about the same amount of energy before reaching a balloon-borne experi¬ 
ment like GAPS. However, the atmospheric production is also important for the 
space-based AMS-02 experiment. Atmospherically produced antideuterons can 
upscatter into space and will create low-energy antideuterons at high geomag¬ 
netic cutoff locations. These effects have been studied semi-analytically in m, 
but atmospheric antideuteron production and transport through the atmosphere 
and geomagnetic field to the ISS should be studied more in the future. 


4 . 5 . Antideuteron interactions in particle physics detectors 

One of the standard tools for studying particle interactions with particle 
physics detectors is the Geant4 simulation suite |19311194j . As mentioned in 
the last two sections, it was also used for the studies of the geomagnetic and 
atmospheric influence. Until recently Geant4 did not allow for the study of light 
antinucleus-nucleus interactions like antideuterons or antihelium with detector 
material. The authors of m added light antinuclei capabilities to Geant4 
using the Glauber approach for the cross sections, the quark-gluon string model 
for annihilations and meson production, and the binary cascade for secondary 
interactions of low-energy mesons. This model was validated between 0.1 GeV jn 
and 1 TeV /n. However, it has to be kept in mind that only very little data on the 
interactions of antideuterons with targets exist. It is key to apply the Geant4 
extended capabilities to geomagnetic, atmospheric, and detector simulations in 
the future and to study the uncertainties in the cross section systematically. 


4 . 6 . Extragalactic transport of antideuterons 

While the focus of the section on propagation was so far on antideuterons 
originating in our Galaxy, one might also wonder if it were possible to search 
for a small component of antideuterons that originate in other galaxies, which 
could then be employed as a test of universal baryon symmetry. The short 
answer is no, intergalactic magnetic fields provide an insurmountable barrier to 
the transport of antideuterons from cosmological distances at typical cosmic-ray 
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energies and limits on diffuse Y-ray emission provide strong constraints on large 
domains of primordial antimatter anywhere within the observable universe. 

The arguments follow those that have been used in the past to argue against 
using cosmic-ray antihelium as a probe of primordial antimatter |1961 1197] . 
Galactic magnetic fields are typically a few p-G. At the intergalactic scale of 
about 1 Mpc, one would naively expect that ejected galactic magnetic fields 
would contribute to intergalactic fields at the 0(1 nG) level. The leakage of 
the observed large dipolar magnetic fields in galaxies would contribute to the 
intergalactic field at least at the 0(1 pG) level. If disordered, these fields would 
prevent diffusion of 0(1 GeV) cosmic rays from even the closest galaxies. The 
only way around this is if cosmic rays were to follow ordered intergalactic fields 
stretched between galaxies. Under these optimized conditions, Ref. [Ml argued 
that cosmic rays would have a diffusion distance of at most 32 Mpc. Additionally 
it was shown that galactic magnetic fields will result in an accessibility problem. 
If cosmic rays have a large accessibility to galaxies, they will enter them and 
be destroyed at each galactic encounter. If extragalactic cosmic rays have a 
small galactic accessibility, they could in principle avoid destruction by galaxies 
en-route, but then they could not get into our galaxy to be observed at any 
significant level. 

Both |196j and [I97j have shown that the observed diffuse y-ray background 
is inconsistent with the existence of significant antimatter domains within the 
Hubble volume. Together, the extragalactic cosmic-ray transport argument and 
the y-ray argument explain why extragalactic antinuclei heavier than hydrogen 
have not been, nor will ever be observed in the cosmic rays. Any antideuterons 
seen in the cosmic rays at typical GeV to TeV energies would have to originate 
within our Galaxy. 


5. Experiments for the detection of cosmic ray antideuterons 


As the search for cosmic-ray antideuterons is a rare-event search, particle 
identification and background rejection are keys. BESS-Polar II, GAPS, and 
AMS-02 exploit different experimental designs to search for antideuterons. Both 
BESS-Polar II and AMS-02 utilize magnetic spectrometers to distinguish an¬ 
tiparticles, while GAPS uses a novel exotic atom capture and decay technique. 
As with all rare signal searches, these different detection techniques are essential 
to validate any possible detection. 


The BESS Antarctic flight program, discussed in Sec. |5.1[ has already pro¬ 
vided the current best upper limit on the antideuteron flux, and ongoing analysis 
of BESS-Polar II data will soon update this limit. AMS-02, installed on the In¬ 
ternational Space Station, has already collected over four years of cosmic-ray 
data, and it will continue collecting data in the years to come. The AMS-02 
collaboration is currently preparing its first antideuteron results, as detailed in 
Sec. |5.2[ GAPS, described in Sec. |5.3[ completed a successful prototype balloon 
flight in June 2012, and is currently proposing for a first Antarctic balloon cam¬ 
paign. In the coming years, AMS-02 and GAPS should provide either the first 
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Figure 15: Cross sectional view of the BESS-Polar II instrument. 


detection of a cosmic-ray antideuteron or upper limits yielding harsh constraints 
on viable dark matter models. 

5.1. BESS-Polar II: The Balloon-borne Experiment with Superconducting Spec¬ 
trometer 

The BESS-Polar program exploits particle tracking in a solenoidal mag¬ 
netic field to identify antimatter. The original BESS-Polar experiment flew 
over Antarctica in late 2004, providing the current best antideuteron flux upper 
limits. The BESS-Polar II experiment collected about 30 days of Antarctic flight 
data from December 2007 to January 2008. No antideuterons were observed, 
and progress towards an improved flux upper limit is reported here. 

5.1.1. BESS-Polar II instrument design 

BESS-Polar II consists of a large solenoidal magnet, filled by inner drift 
chambers (IDC) and a jet-type drift tracking chamber (JET), and surrounded 
by an aerogel Cherenkov counter (ACC) and a time-of-flight system composed 
of scintillation counter hodoscopes. These components are arranged in a coaxial 
cylindrical geometry, providing a large geometric acceptance of 0.23 m^ sr. 

Fig. [T5| shows a cross sectional view of the BESS-Polar II instrument. A thin 
superconducting solenoid provides a uniform field of 0.8 T. Tracking is performed 
by fitting up to 48 hit points in the JET and 4 hit points in the IDC, resulting 
in a magnetic-rigidity resolution of 0.4% at 1 GV and a maximum detectable 
rigidity of 240 GV. The upper and lower scintillator hodoscopes provide time- 
of-flight (TOF) and dE/dx measurements as well as trigger signals. The timing 
resolution of each hodoscope is 120 ps, resulting in a resolution of 2.5%. The 
threshold-type Cherenkov counter, using a silica aerogel radiator with refractive 
index n = 1.03, can reject electron and muon backgrounds by a factor of 12 000. 
The threshold rigidities for antiproton (p) and antideuteron (d) are 3.8 GV and 
7.6 GV, respectively. In addition, a thin scintillator middle-TOF with timing 
resolution of 320 ps is installed between the central tracker and the solenoid, in 
order to detect low-energy particles that cannot penetrate the magnet wall. 
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5.1.2. BESS-Polar II antideuteron identification 

The event selection criteria used to identify antideuterons follow closely that 
used for the antiproton analysis described in [39]. The major difference is the 
background processes that must be considered. In the antiproton analysis, rel¬ 
ativistic electrons, muons, and pions are the main background sources; for an¬ 
tideuterons, the main background are antiprotons. 

A clear detection of an antideuteron over the plentiful antiproton and 
background events requires reliable measurements of particle rigidity, /3, and 
dE/dx. Special care must be taken to reduce tails of distributions, which are 
caused by scattering, interaction, and misidentified tracks in the detector. In the 
study outlined below, samples of protons and deuterons, as well as antiprotons 
with a mask applied to the antideuteron region, were used to optimize selections 
and to determine efficiencies. 

Below are the requirements used in the preliminary BESS-Polar II analysis 
to isolate antideuteron events: 


Clean single track requirement: First, well-reconstructed, non-interacting 
single track events that pass through the fiducial region are selected. Cuts 
are then applied to the of the track fit, the consistency between track 
and TOF hit information, and the number of 0-hits in the IDC to ensure 
the quality of the track and correct timing information. These cuts also 
eliminate hard scattering events. Events with a large number of hits in 
the tracker that are not associated with a track, which may cause mis- 
measurement of the TOF timing and dE/dx, are also rejected. 

dE/dx requirement: Next, bands of dE/dx vs. rigidity space are identi¬ 
fied to distinguish (anti-)protons, (anti-)deuterons, helium, and c/^/tt 
events. Both positive and negative-curvature events are used, since vari¬ 
ous efficiencies for antideuterons can be estimated using positive-curvature 
deuteron events. dE/dx vs. rigidity selections are applied to the upper 
TOF, lower TOF, and JET chambers. The JET dE/dx has better resolu¬ 
tion than the TOF dE/dx, providing good separation of antiprotons from 
antideuterons (and protons from deuterons) up to about 2 GV. 


ACC requirement: The ACC veto, based on Cherenkov counter measure¬ 
ments and designed originally for antiproton identification, are applied to 
reject relativistic e//x/7r samples. In addition, this cut eliminates antipro¬ 
tons and protons with rigidities higher than the threshold R > 3.8 GV and 
reduces the number of events with interactions outside the tracker. 


Mass requirement: Fig. 16 shows the scatter plot of (3 


-1 


vs. 


rigidity for 

events remaining after the dE/dx and ACC requirements. Clear band 
structures are visible, each of which corresponds to particle mass as m = 
TZy^l/P'^ — 1. The antideuteron band width defines a three standard de¬ 
viation detection region. The antideuteron region is initially masked for 
blind analysis. Antiproton events that survived the dE/dx and ACC cuts 
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Figure 16: Left: Plot of /3“^ vs. rigidity for events remaining after the dE/dx and the ACC 
requirements. Right: View of the region indicated by the blue box to the left. The red and 
blue curves show signal and excluded region, respectively. No antideuteron candidate was 
found in the signal region. 


are clearly observed. These remnants are carefully studied, and the selec¬ 
tion criteria are iteratively optimized to remove this tail. 


To select the final antideuteron events, an antiproton excluded region is de¬ 
fined as indicated in Fig. 16 This region is defined such that less than 0.1 event 
is expected to contaminate the antideuteron signal region, assuming that the 
antiproton distribution is Gaussian. No antideuteron candidate is observed 
in the signal region. 


5.1.3. BESS-Polar II current status 

The BESS-Polar II flight was successfully carried out in December 2007 
through January 2008, near solar minimum period. 4.7 • 10® cosmic-ray events 
were collected at altitudes of 34-38 km (average residual air of 5.8g/cm®). Dur¬ 
ing the 24.5 days of observation, all detectors exhibited the expected perfor¬ 
mance, except for the central tracker which showed an instability due to high- 
voltage fluctuation. However, through the development of a time-dependent 
tracker calibration, more than 90% of the data has been successfully calibrated. 
The data has been analyzed for antideuteron signals, as described above, and 
no candidate events were observed. A publication of this result, including a new 
current best upper-limit on the cosmic antideuteron flux, is in preparation. 


5.2. AMS-02: The Alpha Magnetic Spectrometer 

The Alpha Magnetic Spectrometer (AMS-02) is a general purpose high- 
energy particle physics detector, installed on the International Space Station 
(ISS) since 19 May 2011 [36]. A detector prototype, AMS-01, was successfully 
flown aboard U.S. Space Shuttle Discovery in June 1998 |198j . Although its 
mass and volume are optimized to fit the Space Shuttle payload bay, the lay¬ 
out of the AMS-02 detector is typical of larger particle physics experiments. 
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with particle identification performed using the tracking of charged particles 
in a magnetic field, momentum measurement in a time-of-flight system, and 
energy measurements performed in Cherenkov and electromagnetic calorimeter 
systems. 


5.2.1. AMS-02 instrument design 


The AMS-02 instrument [10911200] is shown in Fig. 17 It consists of: a per¬ 
manent magnet {B = 0.14T) and a 9-layer precision silicon tracker, enabling 
the measurement of rigidities up to about 2 TV and a position resolution of 
10 i^m for single charged particles [201] : a time-of-flight (TOF) detector consist¬ 
ing of four detection planes, providing a time resolution of 120 ps for Z = 1; a 
transition radiation detector (TRD) with 20 layers, with the primary goal of e/p 
separation at energies up to hundreds of GeV |202j : a ring imaging Cherenkov 
detector (RICH) with a dual radiator, consisting of silica aerogel {n = 1.05) 
and sodium fluoride (NaF, n = 1.33), for accurate velocity measurement at 
P > 0.95 for aerogel and /3 > 0.75 for NaF [203] : an electromagnetic calorime¬ 
ter (ECAL) composed of nine superlayers with a total thickness of 17 radiation 
lengths, for e'^/e~ and photon measurement 


and anti-coincidence counters 
(ACC) [205] . The absolute value of particle charge is measured independently 
by the tracker, TOF, TRD, and RICH, while its sign is given by the combi¬ 
nation of curvature information from the tracker and directional information 
from the TOF. The acceptance of AMS-02 (without requiring ECAL crossing) 
is approximately 0.5 m^ sr 


5.2.2. AMS-02 antideuteron identification 

A key element in antideuteron identification and background exclusion is 
the determination of particle mass. The mass is calculated from the particle’s 
charge Z, its rigidity TZ and its velocity /?: 


m = Z^Z^ / ^ - 1 • 


(24) 


Three kinetic energy regions can be considered for potential antideuteron iden¬ 
tification. The lowest energy region, Ekm.TOF ~ 0.2-0.6 GeV/n, is the most 
relevant for dark matter studies. In this region, the velocity measurement is 
provided by the TOF. Two higher-energy regions use the RIGH detectors, and 
are thus limited by the Gherenkov thresholds of its radiators: Ekin.NaF ~ 0.6- 
3.8GeV/n for NaF and i?kin,Agi ~ 2.9-8.5GeV/n for aerogel. The NaF radiator 
covers about 10% of the AMS-02 acceptance |207j ]. but, compared to the TOF 
region, this is partly offset by the lower geomagnetic deflection (Fig. 13). 


Nearly all major contributions to the cosmic-ray flux must be considered 
as potential backgrounds to the antideuteron search. The most challenging 
backgrounds are those from particles that have the same mass, charge, or 
mass/charge ratio as the antideuteron, or that may easily produce such a par¬ 
ticle through interaction with the detector. 

AMS-02 data contain many valuable variables for exclusion of these back¬ 
grounds. Selecting clean tracks without in-flight interactions is essential. A 
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Figure 17: Cross sectional view of the AMS-02 experiment. 


large number of hits in the ACC or off the main track is a typical signature of 
a noisy event or pile up of coincident events. Additionally, a match of TRD 
and tracker tracks and a minimum number of TOF hits may be required. In 
the lowest energy region, which uses TOF velocity for mass reconstruction, the 
fact that antideuterons are significantly slower than minimum ionizing particles 
leads to a higher average energy deposition in both ECAL and TOF detectors 
in comparison with lighter particles in the same energy range. 

Specific features of different particle backgrounds are presented below. 

Antiprotons may present the greatest challenge, since they can only be sep¬ 
arated from antideuterons by accurate mass measurement. As mentioned 
above, particle mass is derived from the combination of measured rigidity, 
velocity and charge. The AMS-02 collaboration is currently addressing the 
closely related problem of d/p separation (Fig, [l^, and this analysis will 
serve as a guide for the d/p analysis. A rejection factor of about lO'^-lO® 
will be required. 

Electrons are the most abundant cosmic rays with negative charge. However, 
they can be discriminated from antideuterons using their much lighter 
mass and non-hadronic interaction signatures. Virtually all electrons mea¬ 
sured by AMS-02 have /3 « 1, meaning that only those with a poor veloc¬ 
ity reconstruction (/3 ^ 1) will pose a problem. Additional discrimination 
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Figure 18: Distribution of reconstructed mass in AMS-02 for simulated proton and deuteron 
events for a reconstructed kinetic energy range of 1-2 GeV. In each case, the top distribution 
includes all events with a single reconstructed particle, while subsequent distributions show 
the result of applying quality cuts leaving 60% and 30% of the original well-reconstructed 
deuteron sample (defined as having mrec within ±20% of rrid). 


may be obtained using the energy deposition pattern in the TRD and the 
shower shape and total energy deposition in the ECAL. A rejection factor 
of about 10®-10® will be required. 

Positrons, being less abundant than electrons and with an opposite charge, 
should be excluded efficiently using the same procedures applied to elec¬ 
trons. A rejection factor of about 10®-10^ will be required. 

Protons and deuterons have positive charge, meaning that any backgrounds 
due to these particles must be caused by charge signed misidentification. 
This may happen in AMS-02 if the rigidity is reconstructed in such a 
way that it has the wrong curvature. It could also occur if an upward 
traveling particle is misreconstructed as downward traveling due to an 
incorrect TOF measurement, or if an interaction occurs with a primary 
particle generating an upward traveling secondary. Although deuterons 
have the same mass as antideuterons, it is very unlikely for a deuteron 
that is misreconstructed in rigidity, velocity, or both to be reconstructed 
at the right mass. The masses of events with misreconstructed charge 
are expected to follow a continuous mass distribution rather than a peak 
at the (anti)deuteron mass. A rejection factor of about 10®-10^° will be 
required. 
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Helium is by far the most abundant cosmic-ray component with Z > 1. Helium 
nuclei frequently interact with the AMS-02 detector, producing secondary 
protons and deuterons that may be misidentified as antideuterons. The 
exclusion of such fragmentation events is therefore the main concern for 
this case. The identification of such events typically relies on a search 
for events containing a large number of particles crossing the detector, 
inconsistencies between charge measurements at different points of the 
detector, or an unexpected trajectory change. However, helium-induced 
background should be taken care of when suppressing particles with a 
positive charge (as discussed above). A rejection factor of about lO'^-lO® 
will be required. 

5.2.3. AMS-02 current status 

AMS-02 has been acquiring data on the International Space Station at an 
average rate of approximately 600 Hz [36] . The total number of events detected 
as of 19 February 2015 is over 6.1T0^°, with data collection expected to continue 
for a total period of about 20 years |36j . 

The sensitivity of AMS-02 to antideuterons has been addressed in the past |S2| . 
However, this analysis was conducted in 2007, during the development phase of 
the experiment, and was based on the assumption of a superconducting magnet 
configuration with B — 0.86 T instead of the permanent magnet’s 0.14 T and a 
substantially different tracker plane layout spanning about 1 m instead of about 
2.5 m. The projected antideuteron sensitivity of AMS-02 shown in Fig. is 
based on the superconducting-magnet configuration. An updated study based 
on the final configuration is therefore needed, and will profit from knowledge 
gained in recent years from improved simulations and detector operation. 

Studies of particle identification are currently being performed using both 
AMS-02 data and large Monte Carlo samples produced using a detailed simu¬ 
lation of the AMS-02 detector based on the Geant4 toolkit [19311194] . This is a 
continuous process where ISS data inform the detector construction and physics 
lists used in the simulations. For example, the high statistics of AMS-02 data 
mean that nuclear cross sections may become the dominant source of uncer¬ 
tainty for cosmic-ray predictions, which in turn are assessed using cosmic-ray 
flux ratios such as d/^He [208] . Other ratios like d/'^He may provide further 
constraints on propagation parameters already tuned to B/C data, allowing 
a better estimate of sensitivities for antiproton and antideuteron dark matter 
searches nog. 

Antiprotons will be the most challenging background to antideuteron searches. 
Fig. |18| demonstrates the narrowing of the proton and deuteron mass peaks for 
the reconstructed kinetic energy range of 1-2 GeV, shown after various stages of 
quality cuts (number of ACC hits, track fit quality, number of TOF hits, etc.) 
for simulated data. The selections reduce significantly the proton high mass tail 
while maintaining a high deuteron efficiency. In addition, a second peak visible 
in the deuteron case is due to protons produced in detector interactions. Thus a 
solid understanding of deuteron interactions with the detector will be required. 
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which can be studied with experiments like NA61/SHINE (Sec. 3.2.3). The h- 
nal antideuteron selection cut will ultimately be in a range around the deuteron 
mass peak. In simulated events, a mass resolution iTm/m ~ 12% has been at¬ 
tained for protons at kinetic energies E Ki 0.5 GeV. The AMS-02 deuteron and 
antideuteron analyses are currently ongoing. 

The AMS-02 collaboration has already published results with unprecedented 
accuracy for fluxes and/or ratios on several of the most abundant particle species 


{p [210], He mu, e+/e- 


1^ . B/C |212j l. Future publications are 


expected to include data on antiprotons, deuterons, and heavier nuclei. 


5.3. GAPS: The General Antiparticle Spectrometer 

The General Antiparticle Spectrometer (GAPS) experiment is a balloon- 
borne instrument optimized specihcally for antideuteron and antiproton detec¬ 
tion in an unprecedented low-energy range, providing essential complementar¬ 
ity to the AMS-02 antideuteron search. Traditionally, detectors have relied on 
magnetic spectrometers to discriminate between matter and antimatter parti¬ 
cles. GAPS, however, utilizes a complementary technique that relies on the 
production and decay of short-lived exotic atom [51]. In GAPS, a time-of-flight 
detector measures the momentum of an incident low-energy particle, which then 
loses energy passing through layers of semiconducting Si targets/detectors and 
ultimately stops in the target material. A negatively charged antiparticle will 
be trapped by the nucleus of the Si detectors or A1 frame and form an exotic 
atom in an excited state. The exotic atom will then emit X-rays of characteristic 
energy as the antiparticle quickly de-excites, before annihilating on the target 
nucleus, producing a shower of pions and baryons. 

The X-ray energies, annihilation product multiplicity, dE/dx energy loss, 
and stopping depth provide discrimination between antideuterons and its main 
background, antiprotons. An antiproton with the same velocity will penetrate 
less deeply, produce X-rays of different characteristic energies, and generate 
fewer annihilation products. The combination of these multiple signatures gives 
GAPS enormous rejection power against cosmic-ray protons, electrons, or other 
non-antimatter species. GAPS will also provide the first precision measurement 
of the cosmic antiproton spectrum below 0.25 GeV/n, which will help both con¬ 
strain propagation parameters and probe low-mass dark matter and primordial 
black hole models m- 

The GAPS technique has been validated in antiproton beam tests ina, and 
a prototype balloon instrument (pGAPS) was successfully flown from Taiki, 
Japan in June 2012 [7(11 1192j . Since a magnet is not needed and the GAPS in¬ 
strument is relatively simple, a balloon or satellite payload of significantly larger 
geometric acceptance is possible within a realistic mass budget. In addition, the 
GAPS trajectory will be close to the lowest-possible geomagnetic cutoff, dra¬ 


matically increasing the number of detectable low-energy antiparticles (Fig. 13). 


GAPS is foreseen to execute several Long Duration Balloon (LDB) campaigns 
from Antarctica. 
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Figure 19: Mechanical drawing of the GAPS balloon instrument, with ten layers of Si(Li) 
detectors (gray) surrounded by a plastic scintillator TOF (blue). 


5.3.1. GAPS instrument design 

GAPS consists of a thin plastic scintillator TOF surrounding a 1.6 m x 1.6 m 
X 2 m cube of ten planes of lithium-drifted silicon (Si(Li)) detectors, as shown in 
Fig.p^ The scintillator paddles will be 0.5 cm thick and 160-180 cm long. Each 
paddle will be read out on both ends, providing 500 ps timing resolution. The 
Si(Li) detectors will be 2.5 mm thick with a diameter of either 5 or 10 cm, pro¬ 
viding a total of 10.5 m^ of active area. The Si(Li) detectors serve as degrader; 
exotic atom target; tracker for the incident charged particle and annihilation 
products; and X-ray energy spectrometers. A strip size of 20 cm^ provides suf¬ 
ficient spatial resolution to differentiate between separate X-rays, annihilation 
products, and incident background particles. A readout with both a low-gain 
and high-gain channel will provide both sensitivity to heavy annihilation prod¬ 
ucts and the 4 keV X-ray energy resolution necessary to distinguish antiprotonic 
from antideuteronic exotic atoms. 

GAPS will have a geometric acceptance of about 18 m^ sr in the 0.05-0.25 GeV jn 
energy range. As shown in Fig. the payload will offer a sensitivity of about 
2 • 10“® m“^ s“^ sr“^ (GeV/n)“^ over the above energy range for 105 days of 
flight, achievable with three Antarctic LDB flights of modest duration. With 
this integration time and energy range, the antideuteron measurement is essen¬ 
tially background-free, with about 0.01 expected antiprotons misidentified as 
antideuterons and a cosmic secondary/tertiary rate several orders of magnitude 
below the predicted sensitivity [54]. The GAPS sensitivity to antideuterons is 
an improvement of more than two orders of magnitude over the current limits 
reported by the BESS collaboration. 

5.3.2. GAPS antideuteron identification 

To estimate the sensitivity of GAPS to an antideuteron signal, a two-part 
Monte Garlo simulation based on Geant4 [TMITM] was conducted [54]. The 
first simulation estimated the grasp and flux of stopped antideuterons, including 
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energy loss and in-flight annihilation in the atmosphere and the instrument. 

The second simulation estimated the energy spectrum in the detector due to 
the decay of the exotic atoms, including all instrument interactions. 

The grasp (T) is defined as the product of the geometrical acceptance (G) 
and the stopping efficiency. The grasp was estimated from simulations of an- 
tideuteron interactions as implemented in Geant4. This simulation predicts 
/ Tpdi? = 0.15(0.45) m^ sr GeV/n for antiprotons in the energy range 50-110 
(110-250) MeV/n at the top of the atmosphere, and / Tjdi? = 0.16(0.41) m^ srGeV/n 
for antideuterons in the energy range 50-110 (110-250) MeV/n. 

The second simulation predicts the TOF and dE/dx of the incident an¬ 
tiparticle, as well as the X-ray and nuclear annihilation product energy spectra. 

Since only low-energy antiprotons will produce a signature of an incoming par¬ 
ticle with X-rays and decay products originating from a vertex, they are the 
only background considered here. The antiproton flux is more than four orders 
of magnitude larger than the predicted dark matter antideuteron flux. Thus an 
antiproton rejection factor of more than 10® is required to obtain a 99% confi¬ 
dence level for the detection of two antideuterons. Although, protons are the 
most abundant cosmic-ray species they are not able to fake the distinct exotic 
atom signature of antideuterons and antiprotons. 

The following variables can be used to distinguish antideuterons from an¬ 
tiprotons and is described in detail in [5i] : 

X-ray energies: An antideuteron exotic atom formed in a Si target will pro¬ 
duce X-rays with energies 30keV, 44keV and 67keV; an antiproton will 
produce X-rays with energies 35keV, 58keV, and 107 keV. The X-ray 
yields were estimated with the simple cascade model discussed in |214j . 

The antiproton rejection factor due to one or more antideuteronic atomic 
X-ray detections is about 40, while the antideuteron efficiency is «10%. 

Nuclear annihilation product multiplicity: The intranuclear cascade frame¬ 
work [21511216j can be applied to the annihilation of the antideuteron with 
the Si nucleus, resulting in two different models. The first model assumes 
that the two antinucleons of the antideuteron interact with the Si nucle¬ 
ons simultaneously, while the second model assumes that the antinucleons 
interact with nucleons separately. Requiring a pion multiplicity of five or 
more provides an antiproton rejection factor of about 15 with an average 
antideuteron efficiency of about 60% for the two models. The intranuclear 
cascade model also predicts proton and neutron production. For antipro¬ 
tons incident on Si, the average multiplicity of protons with energy larger 
than 60MeV, chosen to assure good track reconstruction from passage 
through three or more Si(Li) layers, is (Mp) = 0.37. For antideuterons, 
the average multiplicity is (Mp) = 2.06 for the average of the two models. 

A proton multiplicity requirement (Mp) > 3 yields an antiproton rejection 
factor of 150 and an average antideuteron efficiency of about 35% for the 
two models. 

Stopping depth and dE/dx energy loss: The stopping range for antideute- 
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rons is approximately twice as large as that for antiprotons with the same 
TOF. Therefore, antideuterons stopped in the same layer (same stopping 
range) can have a longer time-of-flight. For instance, antiprotons with 
the incoming angle of (0 ± 5)° ((45 ± 5)°) and stopped in layer 3 will be 
rejected by a factor of more than a 1000 with an antideuteron selection 
cut, t'TOF > 13.3 ns (17.6 ns), while the antideuteron efhciency is at about 
80%. Note that the stopped position of the incoming antiparticle can 
be determined by tracking the annihilation products back to a common 
vertex. Moreover, dE/dX energy deposits in the Si(Li) detector can be 
also used as a powerful discriminator for any possible signal event since 
antideuterons at the same velocity as antiprotons have roughly twice as 
much kinetic energy, and deposit more energy in the Si(Li) detector be¬ 
fore stopping. Also the different change of dE/dx from layer-to-layer for 
antideuterons and antiprotons of the same velocity provides additional 
rejection power. In the future, this will be evaluated further with more 
detailed simulations and analyses, such as an event-by-event analysis and 
likelihood analysis. 


5.3.3. GAPS current status 

All major mission challenges for the GAPS design have been solved. Op¬ 
timization of the payload design is now being refined through detailed simula¬ 
tions. Si(Li) detectors with the requisite performance have been produced using 
in-house fabrication techniques. ISAS/JAXA has also demonstrated a cooling 
system using a passive Oscillating Heat Pipe (OHP) technique [2171121811219] . 
which provides the cooling necessary for the Si(Li) system at a reduced mass 
and power requirement compared to standard active pump systems. 

To validate the GAPS instrument concept, a prototype GAPS payload was 
constructed and flown from the Taiki Aerospace Research Field in Japan on 
3 June, 2012 [TOl [TlJ I192| . The payload, shown in Fig. 20 consisted of six 


commercially-acquired lOcm-diameter Si(Li) detectors arranged in three planes, 
with two layers of plastic-scintillator TOF above and one layer below. Because of 
the limited geometric acceptance, stopping power, and flight duration, pGAPS 
did not have any sensitivity to cosmic-ray antiparticles. The goals of the flight 
were instead to verify operation of the Si(Li) detectors in a balloon environment, 
measure the velocity of cosmic rays using the prototype TOF, verify X-ray and 
cosmic-ray backgrounds at flight altitude, and validate both the thermal model 
of the active pump cooling system and the prototype OHP system. 

The total flight duration was about 5 hours, with more than 3 hours at float 
altitude. For about 245 minutes, the payload operated in TOF trigger mode, 
with an additional about 50 minutes devoted to repeated illumination with an 
X-ray calibration tube and a total of about 29 minutes using a Si(Li) trigger 
mode for the study of incoherent X-ray backgrounds. 

All engineering and science goals of the flight were satisfied. Over 600 000 
cosmic-ray triggers were recorded. The thermal model was fully validated, and 
the OHP test was also successful. Both the TOF and the Si(Li) detectors per¬ 
formed in a very stable fashion, with performance as expected according to 
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Figure 20: Photograph of the pGAPS instrument (left) and schematic view of a reconstructed 
cosmic-ray event (right). 


pre-flight measurements. Finally, the incoherent X-ray and cosmic-ray back¬ 
grounds were measured and present no issues to the operation of the full GAPS 
science mission. 

5 . 4 . Experimental path forward 

The use of multiple, complementary experiments has been successfully em¬ 
ployed by the direct dark matter detection community. Separate experimental 
designs yield different backgrounds and approaches to suppressing these back¬ 
grounds, allowing for independent confirmation of any observed signal. As a rare 
event search, the hunt for cosmic-ray antideuterons requires a similar approach. 
If AMS-02 sees one or more antideuteron events a confirming experiment is abso¬ 
lutely needed. GAPS and AMS-02 have complementary energy ranges, but also 
some overlap at low energy, allowing the study of a large energy range for con¬ 
firming signals and the best chance for controlling systematic effects. The key 
virtue comes from the different antideuteron identification techniques. AMS-02 
and BESS rely on magnetic spectrometers, and thus face different backgrounds 
than the GAPS exotic atom approach. A magnetic spectrometer does not only 
have to deliver the isotopic separation between antideuterons and antiprotons, 
but also between the very abundant protons and antideuterons. Limited mag¬ 
netic field strength and tracking resolution can cause protons to be misrecon- 
structed with negative charge and antideuteron mass. In contrast, in the GAPS 
experiment low-energy protons are not able to fake the exotic-atom annihila¬ 
tion signature because protons cannot replace a shell electron nor annihilate 
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with the nucleus. A critical benchmark to reduce systematic uncertainties of 
AMS-02 and GAPS will be the comparison of low-energy (anti)proton fluxes as 
well as antiproton-to-proton ratios. 

AMS-02 orbits at relatively high geomagnetic cutoffs, which reduces the 
number of detectable low-energy charged particles. By flying in Antarctica, 
the proposed GAPS trajectory is specifically tailored to low-energy particles, 
and thus GAPS will face a smaller geomagnetic cutoff correction. In addition, 
AMS-02 was launched at the beginning of the solar activity maximum. Although 
solar cycle 24 exhibits a relatively weak maximum compared to former cycles, a 
first GAPS flight would presumably happen during the next solar activity min¬ 
imum mo], easing the uncertainties associated with low-energy measurements. 

Therefore, the combination of AMS-02 and GAPS antideuteron searches is 
highly desirable. Unless GAPS is begun soon, confirmation of AMS-02 results 
would be a long while forthcoming. Similarly, a non-detection by AMS-02 may 
simply mean that the high background in the AMS-02 orbit and the difficulty 
of rejecting backgrounds are indicative not of the absence of antideuterons, but 
of the presence of high background. 

6. Conclusion 

This article reviewed how antideuterons may be generated in dark mat¬ 
ter annihilations or decays, offering a potential breakthrough in unexplored 
phase space for dark matter. The unique strength of a search for low-energy 
antideuterons lies in the ultra-low astrophysical background for this channel. 
Many dark matter models are capable of producing an antideuteron flux that 
is within the reach of the space-based AMS-02 and balloon-borne GAPS. The 
combination of these two experiments allows for independent experimental con¬ 
firmation, which is critical for a rare event search such as the hunt for cosmic- 
ray antideuterons. After reviewing different classes of models, it is evident that 
cosmic-ray antideuterons are sensitive to a wide range of theoretical models, 
probing dark matter masses from 0(1 GeV) to 0{1 TeV). 

The important uncertainty of the antideuteron production was also exten¬ 
sively discussed. The standard approach is to describe the merging of the an- 
tiproton and antineutron into an antideuteron in coalescence models. These 
coalescence models are tuned to collider measurements of deuteron and an¬ 
tideuteron production. However, it is currently an open question if the an¬ 
tideuteron production depends on the exact underlying process and on the avail¬ 
able center-of-mass energy or if Monte Garlo generators need further refinement. 
The predicted cosmic antideuteron flux arriving at the Earth also relies on the 
modeling of charged-particle transport both in the Galactic medium and in the 
heliosphere. The modeling of Galactic transport is the main source of theoret¬ 
ical uncertainty in the low-energy antideuteron range, while solar modulation 
is more relevant to shaping the low-energy tail of the predicted flux. Galac¬ 
tic propagation is an important uncertainty in predicting antideuteron fluxes 
at Earth. Recent positron data exclude the MIN Galactic propagation model, 
which predicts the lowest antideuteron flux levels at Earth, supporting higher 
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antideuteron flux predictions that could be detected by GAPS or AMS-02. After 
arriving in the solar system, antideuterons can be deflected away from balloon- 
borne or space-based detectors by the geomagnetic field. After traveling through 
the geomagnetic field antideuterons can interact with Earth’s atmosphere. The 
latter two processes can be understood in the same framework and show that 
the choice of geomagnetic location is crucial, but that the influence of the at¬ 
mosphere is a well-controlled systematic effect if a balloon experiment’s flight 
altitude is high enough. 

Cosmic antideuterons have not yet been detected, with the current best flux 
upper limits provided by the BESS experiment. More sensitivity will be pro¬ 
vided by the AMS-02 experiment, which is currently taking data onboard of 
the International Space Station, and by the General Antiparticle Spectrometer 
(GAPS) experiment, which is planned for Antarctic balloon campaigns. The 
combination of AMS-02 and GAPS covers a large energy range with different 
detection methods and backgrounds. Both BESS and AMS-02 utilize magnetic 
spectrometers to distinguish antiparticles, while GAPS uses a novel exotic atom 
capture and decay technique. The BESS collaboration will soon update their 
antideuteron limits. AMS-02 has already collected more than four years of 
cosmic-ray data, and it will continue collecting data until the end of the In¬ 
ternational Space Station. The AMS-02 collaboration is currently preparing its 
first antideuteron results. GAPS completed a successful prototype balloon flight 
in June 2012, and is currently proposing a first Antarctic balloon campaign. In 
the coming years, AMS-02 and GAPS should provide either the first detection 
of a cosmic-ray antideuteron or upper limits yielding harsh constraints on viable 
dark matter models. As with all rare signal searches, complementary experi¬ 
ments will be essential to validate any possible detection. Hence, searching for 
antideuterons with both AMS-02 and GAPS is very important. 

In conclusion, with intensified activities in many related areas cosmic-ray 
antideuterons move away from being the most-unexplored indirect dark matter 
detection technique and are on the way to become an important additional 
channel with breakthrough potential for disentangling the nature of dark matter. 
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